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Abstract 

We consider the question of conformal invariance of the long-range Ising model 
at the critical point. The continuum description is given in terms of a nonlocal field 
theory, and the absence of a stress tensor invalidates all of the standard arguments 
for the enhancement of scale invariance to conformal invariance. We however show 
that several correlation functions, computed to second order in the epsilon expansion, 
are nontrivially consistent with conformal invariance. We proceed to give a proof of 
conformal invariance to all orders in the epsilon expansion, based on the description 
of the long-range Ising model as a defect theory in an auxiliary higher-dimensional 
space. A detailed review of conformal invariance in the d-dimensional short-range 
Ising model is also included and may be of independent interest. 


August 2015 



Contents 


1 Introduction 3 

2 Field-theoretical setup 5 

3 Tests of conformal invariance 8 

3.1 (003) 9 

3.1.1 Explicit computation. 10 

3.1.2 General argument. 12 

3.2 (0204). 15 

4 Conformal invariance of the SRI critical point 16 

4.1 Scale and conformal invariance breaking effects . 19 

4.1.1 A proof which is too quick, and a paradox. 19 

4.1.2 Correct proof. 20 

4.2 SRI conformal invariance beyond perturbation theory. 23 

5 Conformal invariance of the gaussian phase 25 

5.1 Direct argument. 25 

5.2 Argument from correlation functions. 26 

5.3 AdS argument. 26 

5.4 Caffarelli-Silvestre trick. 27 

6 Conformal invariance of the LRI critical point 29 

6.1 Ward identities. 29 

6.2 From the Ward identities to conformal invariance. 32 

6.3 Bulk-to-defect OPE and the anomalous dimensions. 32 

6.4 Beyond perturbation theory. 34 

7 Conclusions 35 

A The hard diagram 37 

A.l Mellin-Barnes method . 38 

A.2 Gegenbauer polynomial method. 40 

A.2.1 Integration rules . 41 

A.2.2 Computation. 42 

B Conformal invariance of composite operators 44 

C Relative normalization of the 0 and 0^ OPE coefficients 46 


2 






















1 Introduction 


We will be studying the long-range Ising (LRI) model—a close cousin of the usual, short- 
range, Ising model (SRI). While in the SRI only nearest-neighbor spins interact, the LRI 
energy involves interactions at arbitrary distances with a powerlaw potential: 

= ( 1 - 1 ) 

where d is the space dimensionality, and cr > 0 is a free parameter. The sum is over all pairs 
of lattice sites and is the distance between them. We are considering the ferromagnetic 
interaction case J > 0. 



Figure 1: The three phases for the LRI critical point. The boundaries are a = dj2 (straight 
solid) and a = 2 — r/sRi(d) (curved dashed, interpolated using the known exact r]sRi{d) for 
d= 1,2,4 and numerical values from the e-expansion [1] and the conformal bootstrap [2] for 
a few intermediate d). Here we will be working near the boundary a = d/2. 


Just like the SRI, the LRI has a second-order phase transition at a critical temperature 
T = Tc- The critical theory is universal, i.e. independent of the short-distance details such 
as the choice of the lattice. It has however an interesting dependence on a (see Fig. 1): 

1. For a < d/2, the critical point is a gaussian theory described by the nonlocal action 
involving the “fractional Laplacian” operator 

So = j d'^xd'^y(j){x)(j){y)/\x-y\'^^'^ oc j d'^x (j) C^cj). (1.2) 

The field 0 represents the spin density, and its scaling dimension is read off as 

A^ = {d-a)/2. (1.3) 

Composite operators can be bnilt ont of 0 by differentiations and taking normal-ordered 
products. Since the theory is gaussian, there are no anomalous dimensions. 

2. For d/2 < a < a^, the critical point is a nontrivial, non-gaussian, theory, whose held- 
theoretic description can be obtained by perturbing the nonlocal action Sq by a local qnartic 
interaction J d'^x (/"^{x). In the range a > d/2 this interaction is relevant and generates a 
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renormalization group (RG) flow, reaching a fixed point in the IR. This IR hxed point is 
believed to be in the same universality class as the critical point of the LRI lattice model. 
Interestingly, as it will be explained below, the dimension of 0 at the hxed point is still 
given by the same formula (1.3). However, composite operators do get nontrivial anomalous 
dimensions. 

3. Finally, for a > a* the potential is so strongly peaked at short distances that the LRI 
critical point is identical with the SRI critical point, i.e. is independent of a. The value of 
(T* is determined so that the dimension of 0 is continuous across this transition: 

a, = d-2Af^ = 2-T]sRi, (1.4) 

using the usual dehnition of the r] critical exponent. 

The above picture is considered standard since the foundational work of [3, 4]; it has 
been supported by theoretical studies [5, 6] and by Monte-Carlo simulations [7]. 

There is some ongoing debate about what precisely happens near cr = cr*. Some recent 
studies [8, 9] claim observing deviations from the standard picture in this region, while 
others provide further evidence in its support [10, 11].^ In this work we will be focussing on 
the region near the other phase boundary a = d/2, i.e. away from a = cr*. We will therefore 
not try to weigh in on this debate, except for a small comment in the discussion section. 

Problem of conformal invariance 

As we will review in section 4, the SRI critical point enjoys the property of conformal 
invariance. The utility of this symmetry in 2d is long known, as it allows for an exact 
solution of the critical theory [12] via a method called the conformal bootstrap.^ In 3d, an 
exact solution is not yet known, but the conformal bootstrap has recently been used to get 
the world’s most precise numerical determinations of the SRI critical exponents [13].^ 

What about the LRI critical point? In region 3 it’s identical to the SRI so it’s conformal. 
In region 1 it’s described by the nonlocal gaussian theory, whose conformal invariance is 
well known and will be reviewed in section 5. 

Our main goal here will be to show that the LRI critical point is conformal also in the 
non-gaussian region 2. This is harder to ascertain, and as far as we know, this issue has not 
been previously discussed. Here we will present a proof valid to all orders in perturbation 
theory. Having conformal symmetry also in this region is very interesting and useful, paving 
the way for the conformal bootstrap methods.^ Until now, the LRI critical point was studied 

^Some of this recent literature [10, 11] is phrased in terms of an idea of an “effective dimension” Des{o') 
such that the LRI critical point in d dimensions is supposed to be equivalent to the SRI critical point in Dgg 
dimensions. We find this idea highly implausible. The relations among a handful of LRI and SRI critical 
exponents, which led the authors of [10, 11] to postulate this notion, are most likely approximate (see also 
[9]). In our work the idea of the effective dimension is not going to play any role. 

^In fact, the 2d SRI critical point is invariant under the Virasoro algebra, which is an infinite-dimensional 
extension of the global conformal algebra. In this paper, we will be working in general d and by conformal 
invariance we will mean global, finite dimensional, conformal invariance, unless specified otherwise. In 2d 
this is sometimes referred to as Mdbius invariance. 

^For prior work see [14-19]. An alternative technique has been developed in [20-22]. 

^Some work has already been done in [16], section 5.4, and [23]. 
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using RG methods, both perturbative [3-6] and nonperturbative [11]. 

The remainder of this paper is organized as follows. In the next section we discuss 
the basic setup of the epsilon expansion for the long-range Ising model. We then provide 
nontrivial evidence for conformal invariance by computing (00^) and up to order e^, 

as we hnd that these correlators vanish at the hxed point. This behavior does not follow 
from scale invariance alone but it is a necessary condition for conformal invariance, so we 
are led to believe that the LRI at criticality could actually be conformally invariant. The 
remainder of the paper is dedicated to an all-orders proof of this claim. In sections 4 and 
5 we review the proofs of conformal invariance of the SRI and of the gaussian theory (see 
below for a dehnition). This sets the stage for the proof of conformal invariance of the LRI, 
which we present in section 6. In that section we also discuss the prospects for proving 
conformal invariance at the nonperturbative level. We end the paper with some concluding 
remarks, and several technicalities have been relegated to two appendices. 


2 Field-theoretical setup 

In this work we will study the LRI critical point for 1 ^ d < 4, inside the non-gaussian 
region 2, close to the boundary separating it from the gaussian region 1, i.e. for 

a={d + e)/2, 0 < e < 1 (2.1) 

The UV dimension of the 0 held is then 

A<^ = (d-e)/4, (2.2) 

so that the quartic term 0^ is a weakly relevant perturbation. The IR hxed point is 
then accessible in perturbation theory. As is standard [5, 6], we will setup a perturbative 
expansion using the analytic regularization scheme, where Feynman diagrams are considered 
analytic functions of e. This is convenient since it allows one to evaluate integrals without 
introducing an explicit UV cutoh. 

As mentioned in the introduction, the relevant action is 

^ = ^o + ^/d"a:0^ S'o = ^yd‘'a:0£<,0. (2.3) 

The coefficient whose precise value is unimportant, will be hxed so that the free 0 two 
point function is normalized to one: 

(0(x)0(O))5o = . (2.4) 

Hid is close to 4, we should assume a stronger condition e ^ 4 —d, so that we stay closer to the boundary 
between the regions 1 and 2 than to that between 2 and 3. In the opposite case e > 4 — d the structure 
of perturbation theory is modified due to the presence of a weakly irrelevant operator 5^0. See [5, 6] for a 
discussion. 
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Notice that the analytic regularization scheme is mass independent. So we can avoid 
introducing the mass term into the action (2.3). In other words, in this scheme the 

flow which leads to the flxed point has = 0 all along the flow. 

As usual in the held theoretical studies of critical phenomena, it will be extremely useful 
to view the theory (2.3) as regulating the theory with e = 0, where the interaction is 
marginal. As e —)■ 0, computations in theory (2.3) give poles in e. These poles can be 
removed order by order in perturbation theory by multiplicatively renormalizing the terms 
in the action.® 

The gaussian term in (2.3) is not renormalized, because it’s nonlocal, while the diver¬ 
gences are local. In other words, the theory under consideration does not have wavefunction 
renormalization, and the bare and the renormalized held 0 coincide. In particular, the 
anomalous dimension of 0 vanishes. This implies that, to all orders in perturbation theory, 
the dimension of 0 at the IR flxed point is equal to its UV dimension (2.2) [3-5]. 

On the other hand, the coupling constant does require renormalization. The bare 
coupling qq is expressed in terms of the finite, renormalized, coupling g as: 

9o = Z{g,e)gg.\ (2.5) 

where /i is the renormalization scale and Z is the renormalization factor which starts with 
1 and contains an ascending series of poles in e: 

Z(j,£) = l + 53“y^(9)£-‘. (2.6) 

The coefficient fk{g) of the pole is a power series in g starting from 0{g^)] it gets 
contributions from all loop orders larger than or equal to k. 

The 0-function of the theory can be expressed in terms of the single pole coefficient: 

/3{g)=dg/d\ogfx = -eg + g‘^f[{g). (2.7) 

The story is exceedingly similar to the e-expansion for the Wilson-Fisher (WF) fixed point 
(see e.g. the books [28, 29], or [30] for a concise treatment), with an extra simplification 
that there is no wavefunction renormalization. 

Let us compute the 0-function at the first nontrivial, one-loop, order. We need to 
determine the one-loop counterterm to the coupling. Consider the four point function of 0, 
which up to the second order in g is given by the sum of diagrams: 


+ X X + perms 


^Multiplicative renormalizability is most familiar in local quantum field theory, but it holds also for 
nonlocal theories of the kind we are considering here, with the standard proof [24]. For the validity of 
Weinberg’s convergence theorem [25] (that if all subdiagrams are superficially convergent, then the diagram 
is convergent) in this context see the discussion in [26], p.600. For the structure of divergences in theories 
with propagators involving arbitrary powers of (studied in the context of analytic regularization), see 
[27] and the discussion in [5]. 
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The 0{g^) contribution will contain a 1/e pole coming from the region where the two 0^^ 
insertions are near each other. The contribution of this region can be easily extracted using 
the operator product expansion (OPE) of the gaussian UV theory: 

X 0^(0) D (72/|a:|'^-") 0^(0). (2.8) 


This OPE implies that the expansion of the functional integral to the second order will 
contain a short-distance contribution corresponding to the insertion of 0^ times 


72 X - 


X 




dd 


X 


— 


9 ^ S, 






(2,9) 


where /T{d/2) is the volume of the unit sphere in d dimensions. The 1/e pole is 

canceled by adding a coupling counterterm 


6g = {36g^/4\){Sd/e). 


( 2 . 10 ) 


Therefore the expansion of the single pole function fi{g) in (2.6) starts with 

h{g) = Kg + 0{g^), iP = 3S,/2. (2.11) 

The one-loop /3-function is given by 

m = -^9 + Kg^ + 0{g^). ( 2 . 12 ) 

It has a zero at 

g = g, = e/K + Oie‘^). (2.13) 

This zero corresponds to the LRI critical point that we want to study. ' 

Below we will be interested in the correlation functions of the composite operators 0"". 
As usual, we will dehne the renormalized operators [0""] = [0”]g,/i which remain hnite in the 
limit e —)■ 0. They are related to the bare operators 0” by rescaling factors subtracting poles 
in e: 

0" = Z„((/,e)r]. (2.14) 

As indicated, the operators [0”] depend on the scale fi and on the value of the coupling g at 
this scale. On the other hand correlation functions of the bare operators don’t depend on /i 
and g separately, but only on their combination g^. As usual, this gives a CS equation for 
the correlators of the renormalized operators. The anomalous dimension of [0”] is given by 

7n(5') = Z0^02'„/01og/i|go=const • (2.15) 

Notice that since the dimension of 0 is not integer, in the considered theory the 0” operators 
don’t mix with the operators where some 0’s are replaced by derivatives.® 

^ In d = 1, 2, 3 and for sufficiently small e > 0 the existence of this fixed point has been shown rigorously 
[31-34]. See also [35] for a series of rigorous results about the LRI phase transition. We are grateful to 
Abdelmalek Abdesselam and Pronob Mitter for communications concerning these works. 

®This is true for generic d, while for integer d there can be some mixing for high n. E.g. « 1/2 in 
d = 2, and one can trade 4 0’s for 2 derivatives. This could become important for operators starting from 
0®. In this paper we will mostly work with operators up to 0"^ and we will ignore this. 
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As already mentioned above, this theory has no wavefunction renormalization. Therefore 
0 is a finite operator without any rescaling. So Zi = 1 and 71 = 0. The UV dimension of 0 
given in (2.2) will also be its IR dimension.® 

Let’s compute the anomalous dimensions of 0"', n ^ 2, at the lowest nontrivial order. 
The 1/e contribution to the correlation functions of 0"' will come from the nearby 
insertions. In this region we can use the OPE generalizing (2.8): 

0"(O) X (t)‘^{x) D [ 6 n(n — l)/|x|‘’*“'’](/"(0). (2-16) 

Integration in x gives a 1/e pole which must be canceled by rescaling the operator. We get 
Zn = l-Kng/e + , iPn = n(n - l)Sd/4 ^ --iriig) = Kng + 0{g^). (2.17) 

At the IR hxed point we have 

7 n ~ Kng* = [n{n - l)/ 6 ]e + O(e^). (2.18) 

These are the same anomalous dimensions as in the usual e-expansion for the WF hxed 
point (see e.g. [36]). Indeed, the answer at this order is controlled by the combinatorics of 
the Wick contractions in the gaussian UV hxed point, which is the same in the local free 
scalar theory and in the nonlocal theory dehned by our Sq. Of course, it’s not true that all 
computations are the same between the two theories. We will see some examples below. 


3 Tests of conformal invariance 

How can we check if a certain model is conformally invariant or just scale invariant? 
The currently available data on the LRI amount to the anomalous dimensions of scaling 
operators, which are determined with RG methods and also measured on the lattice from 
two point functions. These data cannot distinguish between scale and conformal invariance, 
because both predict the same functional form for the critical two point function: 

{0{x)0{<S)) (x\x\-‘^^° . (3.1) 

Here (P is a generic scalar operator of dimension Aq. 

One celebrated prediction of conformal invariance is the form of a three point function 
[37], but this is harder to compute and to measure on the lattice than the two point function. 
We are not aware of any three point function data for the LRI critical point. 

An easier discriminating variable is the two point function of two different operators. 
Scale invariance predicts that 


(C>i(a;)C>2(0)) = ci2|x 


— Ai—Ao 


(3,2) 


®This fact has also been established non-perturbatively for sufficiently small e > 0 for which the existence 
of the fixed point is rigorously known [33, 34], see note 7. 



while conformal invariance implies [37] the stronger constraint that C 12 = 0 unless Ai = A 2 . 
To be precise, this conclusion is reached if both operators are so-called primary operators, 
i.e. if they are not derivatives of other operators. If the Oi are both of the form (5^)"*(9 for 
the same operator (9, they will of course have a nonzero two point function. This trivial 
case is easy to monitor since the scaling dimensions are different by an even integer. 

To summarize, conformal invariance implies that any two scalar operators whose dimen¬ 
sions are not different by an even integer must have a zero two point function, while scale 
invariance would allow such two point functions to be nonzero. 

The goal of this section will be to study the following two point functions of different 
scalars at the critical point of the LRI: 

(0(a:)0^(O)) and (0^(x)0"^(O)). (3.3) 

Notice that we chose pairs of operators with the same parity under 0 —)■ —(f). If the IR hxed 
point is only scale invariant but not conformal, these correlators could be nonzero. 

We will consider these correlators in the perturbative setup of the previous section. As 
we will see, a nontrivial check requires to go to at least the second order in the coupling 
constant gf, or in the parameter e parametrizing the deviation from marginality. At O(e^), 
scale invariance alone then allows for both of the above correlators to be nonzero. 

As a matter of fact, it turns out that the hrst correlator is a bit special, as it involves 
(f whose anomalous dimension is identically zero, as well as 0 ^ which is related to (f via a 
“nonlocal equation of motion” (nonlocal EOM).^° Using these facts, we will give an all-order 
argument that the correlator {(f(f^) vanishes at the IR hxed point. We will buttress the 
argument by an explicit second-order computation showing that the correlator vanishes at 
0 ( 6 ^). 

We do not know any analogous general argument for the correlator (0^0^). This seems 
to be a truly generic correlator involving helds with unrelated anomalous dimensions. We 
will perform an explicit second-order computation for this correlator, hnding that it also 
vanishes at O(e^). 

These computations lead us to believe that the LRI critical point is in fact conformally 
invariant, and to look for a general proof of this fact. A proof that we found will be discussed 
in the rest of the paper. 


3.1 (0 ./.=’> 

By dimensional considerations, this correlator should behave in the IR as We 

recall that at the WF hxed point C 7 ^ 0 because (f^ is not a primary but is related to (f 
by the equations of motion: 0^ oc d'^cf. In the conformal held theory (CFT) language, (f^ 
is a descendant of the primary 0. At the LRI hxed point, we know for sure that 0^ is not 
a descendant of 0 since the dimensions don’t match. So C* 7 ^ 0 would disprove conformal 
invariance, while (7 = 0 would be evidence in its favor. 

thank Riccardo Rattazzi who emphasized this to us. 
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We will first show that C = 0{e^) by an explicit computation, and next present an 
argument that C* = 0 to all orders in perturbation theory. 


3.1.1 Explicit computation 


We will now present an explicit computation of the correlation function (00^) to the second 
order in the coupling. More precisely, we will consider the correlator 

F{x,g,fx) = (0(x)[0^](O)), (3.4) 

where [ 0 ^] = [ 0 ^]^,/^ is the 0 ^ rescaled by subtracting poles in e as discussed in section 2}^ 
This correlator allows an expansion in powers of the renormalized coupling g without poles 
in e. By dimensional reasons, it should have the form 

F = (3.5) 

where A® = 3 A 0 is the UV dimension of 0^ and / is a function of dimensionless variables 
g and s = fi\x\. 

We will have a Callan-Symanzik (CS) equation expressing the invariance of the theory 
under simultaneous changes of g and fi leaving go hxed: 

+ I3i9)9g + 'y3ig)]F{x, g,g.) = 0. (3.6) 

The only difference from the usual 0^ theory is that the anomalous dimension of 0 does not 
appear since it’s identically zero. Substituting (3.5) into (3.6), we get an equation for /: 

[sds + (3{g)dg + 73 (^)]/(s, ^) = 0. (3.7) 

At large distances (3{g) —)• 0 and 73 ( 5 ') —)■ 73 ( 5 '*), and Eq. (3.7) predicts 

f{s,g) (3.8) 


We thus rederived the result that at large distances the considered correlation function 
should go as where A<^3 = A^g^ + 73 ( 5 '*)- This argument does not determine the 

value of C, which in particular may still turn out to be zero. 

To determine the prefactor C, we have to match the CS equation to hxed-order pertur¬ 
bation theory. Recall that the full solution of Eq. (3.7) can be written in the form: 


f{s,g) = /(^(s,^)) exp 


dlog s'73 (s', ^)) 


(3.9) 


^^Note that since we don’t take the limit e —> 0, operator renormalization is not strictly necessary. The 
rescaling factors Z„{g,e) are finite at finite e, and the 1/e terms in them never even become large, since 
g = 0{e) all along the RG flow. The role of renormalization is rather that of convenience, as it allows to 
cleanly separate the 0(1) effects coming from the powers oi g/e from the effects suppressed by powers of g 
which are truly higher order. This makes the weakly coupled nature of the theory manifest. 
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where g is the “running coupling” solving the following differential equation with a boundary 
condition at s = 1: 

sdsg = -l3{g), g\s=i = g ■ (3.10) 

In particular, at s = 1 the function f{s, g) reduces to f{g)- Once this latter function is hxed 
from perturbation theory, the prefactor C in (3.8) is found as 

C = f{9*)- (3.11) 

We will now show that this vanishes up to order e^. 

For this we would like to extract / up to the second order in g. Consider hrst the 
nonrenormalized correlator 

Fo(a;) = (0(a;)03(o)), (3.12) 

which has no tree-level contribution. Up to the second order in the coupling it is given by 
the sum of two position-space diagramsd^ 



(3.13) 


Both these diagrams are easily evaluated using the following basic integral (which in turn 
can be derived by going to momentum space): 


d^y 


Wa Wb 


r(cbii) 

\x-y\^\y\^ WA+B-d\x\^+^-^' ^ ^ ^ r(kl/2) 


Using this result, we obtain 

Fo{x) = Rigo/\x\‘^~‘^^ + R 2 gl/\x\‘^~^^ + 0{gl ), 

where 


'll) d — e 'Wo d—e 

Ri = -^^ ^ -evr 


Wd-2e 


,Wd-eW3. 5 


,/, r(-d/4)r(d/2) 
r(3d/4) 


^2 = X 




_ ±d—-e 

2 2^ 2^ 


2 Wd-2e'Wd-3e 


97r' 


r(3<i/4) 


( 3 . 14 ) 


( 3 . 15 ) 


( 3 . 16 ) 


The given approximate expressions are the leading ones in the small e limit. Notice that 
Ri = 0(e). This has a simple explanation: in the limit e —)■ 0 the interaction becomes exactly 
marginal, the integral dehning Ri becomes conformal, and it should give zero answer for 
a correlator of two helds of different scaling dimensions by the usual arguments based on 
conformal symmetry. To get a nontrivial check of conformal invariance, we have to go to 
the second order in e, hence to the second order in g, which is what we are doing. 

To get at the function /, we need to replace the coupling go by g via (2.5), (2.11): 


go = Zgn^ = [g + Kg^e ^ + 0{g^)]g.% (3.17) 

^^The signs and combinatorial factors are left implicit. 
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and to use the relation between 0^ and [0^], see (2.14), (2.17). This gives the following 
expression for F to the second order in g\ 

F{x) = {gRis^ + g'^[Ri{K + K^)e-h^ + R 2 S^^]] /\x\'^-\ (3.18) 

We see the structure is in agreement with (3.5). To extract / we set s = 1 and obtaind^ 

f{g) = gR^ + g^[R^{K + K^)e-^ + R^] + 0{g ^). (3.19) 

Since Ri = 0(e) this expression is free of poles in e: as we mentioned above the correlator 
and in particular / should have a regular expansion in g without such poles. Using (3.19) 
and (3.9) we can compute the considered correlator at all distances with e^ accuracy. 

Now using the values of the various constants appearing in (3.19), it’s easy to see that 
to the order that we computed it can be rewritten as 

f{g) = -R^e-\-eg + id/) = -Rie-^(5{g ), (3.20) 

This form of writing makes it manifest that f{g) vanishes at the IR hxed point. 


3.1.2 General argument 

We will now give a general argument that C* = 0 to all orders in perturbation theory. The 
idea is to use the nonlocal EOM of the LRI held theory (2.3)d^ 

MaC^+^ct)^ = Q. (3.21) 

This can be used to express the correlators of 0^ in terms of those of 0. In particular: 

(0^(a;)0(O)) oc C(0(a;)0(O)). (3.22) 

This equation is subtle to use, because the fractional Laplacian is a nonlocal operator. We 
will therefore proceed cautiously. 

First of all we have to understand in some detail the correlator (0(x)0(O)). Since 0 does 
not acquire an anomalous dimension, its two point function has the form 

(0(a;)0(O)) = p{x)\x \-^^‘^, (3.23) 

where 

pi.x) = p{g[s,g)) . (3.24) 

^■^We have analyzed also the leading log s terms in the expansion of G around s = 1 and checked that 
they are consistent with what the solution (3.9) to the CS equation predicts. 

^^Notice that it would be a non-permissible stretch of terminology to call (jP a “nonlocal descendant” of 
0 on the grounds of this equation. Descendants are defined as local derivatives of primaries, and we are not 
aware of any useful generalization of the descendant concept to nonlocal relations. 
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is a function of the running coupling g which we introduced in the previous section. The 
function p can be determined by matching to perturbation theory. We will only need to 
know its rough structural properties. 

Through 0{g‘^) we have two diagrams: 


+ 



(3.25) 


This implies that 


p(g) = 1 + Qg2 + ..., Q = (vr76)r (-f)/r (f) . (3.26) 


We see that the function p{x) approaches a constant at short and long distances: 


p{x) 



X ^ 0, 

1 + O(e^), X ^ oo . 


(3.27) 


In what follows we will also need the asymptotics of the approach to the long distance 
limit. From now on we will £x /i to be a scale pc at which the coupling g ^ g^/2 is roughly 
halfway between the UV and the IR hxed points. The long distance asymptotic behavior 
of the running coupling g is given by: 


g = g.- 0(e5-^'("‘)) (s > 1), (3.28) 

where (5\g^) = e + O(e^). It follows that p{x) approaches the long distance limit as 

p{x) = p{g^) + . (3.29) 

We are now ready to compute the long distance behavior of the correlator (0^(x)0(O)). 
According to Eq. (3.22) we need to evaluate the integral: 

/ = j d^yT{x-y){ct>{y)m). (3-30) 

where T{x — y) oc |x — y\~‘^~'^ is the position space kernel of the fractional Laplacian, see 
eq. (1.2). Notice that this kernel is not absolutely integrable and at short distances it must 
be understood in the sense of distributions, as the Fourier transform of \k\^- 

We split the above integral into two parts / = Ii + I 2 , one against the long distance 
asymptotics of and the rest: 

h = j d%T{x - y)p{g.)\y \-^^^, (3.31) 

h = I-h = j d^ynx-y)[p{y)-p{g,)]\y\-^^K (3.32) 
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The first part Ji vanishes at non-coincident points, Ji oc S{x), by the dehnition of Green’s 
function of the gaussian theory. As to I 2 , in the limit of very large x the leading asymptotics 
of the integral will come from large y, where we can use the asymptotics (3.29). We get: 

/2 ~ y d'^yT{x — . (3.33) 

By dimensional analysis, this integral behaves oc |x|“" with 

O' = (T 2A(^ -f- . (3.34) 

It’s important that this exponent is larger than -|- A^s: 

~ (^0 + ^03) = , (3.35) 

using the known anomalous dimension of (j)^, 'jsig*) = € + O(e^). This means that the 
constant C in the natural dimensional asymptotics of the studied correlator, see (3.8), has 
to vanish. The current argument establishes this fact to all orders in perturbation theory. 

The above reasoning was made possible by two properties: the nonlocal EOM and the 
vanishing anomalous dimension of 0. Before leaving this section, we will similarly prove one 
more interesting fact: that 73 = e to all orders in perturbation theory. In other words the 
leading order result in (2.18) is in fact exact for n = 3. 

We will use the nonlocal EOM to prove the following relation between the IR dimensions 
of 0 and 0 ^: 

A^s = A^ + a (IR), (3.36) 

of which 73 = e is an immediate consequence. Eq. (3.36) looks similar to the relation 
A^3 = A^ -|- 2 valid at the WE fixed point, also a consequence of the corresponding EOM. 
However, due to nonlocality, the proof is a bit more subtle for the LRI. 

The idea is to use the nonlocal EOM twice, expressing (0^0^) in terms of (00). This 
relation takes the form: 

(03^3^ = 36^0W^£<,£<,(00)pert , (00)pert = (00) “ (00)so • (3-37) 

with the two acting on each of the arguments of (00)pert- The subtlety here is that we 
have to subtract the gaussian two point function. This is easy to understand in perturbation 
theory. ( 00 ) pert is the sum of all diagrams in which each 0 is connected to a vertex, like in 
the second diagram in (3.25). When we act with the legs connecting 0’s to the vertices 
get cancelled, and we reproduce all diagrams for (0303). Were we to keep ( 00 ) 5 ^, we would 
get an extra nonlocal contribution which does not correspond to any diagram. This is in 
contrast to what happens when using the equation of motion in the local theory. In that 
case the contribution from the unperturbed propagator is zero at noncoincident points and 
we don’t have to worry about it.^^ 

^®One can also check that (3.37) has a smooth limit when go 0, reducing to {(l)^(j)^)so in this limit. This 
would not be the case were we to keep (00) s^. 
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Let us proceed now to the proof of (3.36). It will be convenient to work in momentum 
space. The asymptotics (3.29) of the function p{x) means that the Fourier transform of 
(00)pert behaves at small momenta as: 

(0(-fc)0(fc))pert ~ \kr[l + 0(|fc|^'(^*))]. (3.38) 

It’s important that the proportionality coefficient here is nonzero (we computed that it’s 
0(e^)). When we act on this correlator once by the fractional Laplacian, multiplying by 
the leading term gives 1, which is a delta-function in the position space, and the subleading 
term determines the long-distance asymptotics. We thus reproduce the above result that 
{(j)^ 0) vanishes in the IR faster than its natural scaling. However, when we act by the 
fractional Laplacian twice, the leading result is nonanalytic and it is this result which 
determines the long-distance asymptotics. We conclude that the two point function (0^ 0^) 
behaves at long distances as which can be expressed as Eq. (3.36) . 


3.2 ( 0204 ) 


In the previous section we studied the correlator (00^) and found that it vanishes at the 
IR fixed point, consistently with conformal invariance. However, as we saw, the studied 
correlator was actually a bit special, since it involved two fields related by the nonlocal 
EOM of the LRI. One could wonder if the vanishing of this correlator is an accident. In 
this section we will study the correlator (0^ 0^), which as far as we can see is a truly 
generic correlator of our theory. We will find, by an explicit computation at O(e^), that this 
correlator also vanishes at the IR fixed point. We consider this a strong piece of evidence 
for the conformal invariance of the LRI critical point. 

The computation proceeds similarly to (00^), so we will be brief. The renormalized 
correlator 

H{x,g,p) = ([0^](x) [0^](O)) = h(s,^,/r)|a;rV“V (3.39) 

satisfies a CS equation which can be solved to give: 


his,g,p) = h{g{s,g))expl- f d\ogs' [ 72 (^(s', ^)) + 74 (^(s', ^))] | • (3.40) 


The function h is determined by matching to perturbation theory. We have one diagram at 
0{go) and three diagrams at 0{gl)-. 


+ 





( 3 . 41 ) 


The first two diagrams here are the same as in (3.13) times an extra propagator. The third 
diagram is new, but it can also be readily evaluated using Eq. (3.14) repeatedly. The final 
diagram is the hardest—it is analyzed in appendix A. 
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The nonrenormalized correlator is thus given by 


Ho{x) 


Pl9o {P2a + P2b + P2c)9o , n(n^\ 

^|3(i/2-5e/2 |^|3d/2-7e/2 '-^VVoJ ) 


(3.42) 


where we split the coefficients diagram by diagram. Taking into account signs and symmetry 
factors, we have (see (3.16)): 


Pi — 8i?i, 

P 2 a = 8 R 2 , 


P26 = 4 


Wd-2eW3^_7^ 


Svr^ 


1) 


r(¥) ’ 


P2.«24ir‘‘r (-3/r(f). 


( 3 . 43 ) 


To pass from this result to the function h, we need to perform the coupling and operator 
renormalization. We get: 

^{ 9 ) = 9 P 1 + 9^[Pi{P + P 2 + -^ 4 )^ ^ + P2a + P2b + -P2c] + 0{g^). (3.44) 


Putting in the values of all coefficients, we hnd that this expression is proportional to the 
one-loop beta function, similarly to (3.20). We conclude that h{g) vanishes at the IR fixed 
point at O(e^). 


4 Conformal invariance of the SRI critical point 


In the previous section we found that two point functions (0 0^) and (0^ 0^) vanish at the 
critical point of the LRI model. Nonzero values of these correlators would be consistent with 
the scale invariance and the Z 2 symmetry of the IR hxed point. We know of no symmetry 
but conformal invariance which could forbid these two correlators. We are led to believe^® 
that the model is conformally invariant, and we must understand why. 

The phase diagram of the LRI model in Fig. 1 has three regions. Here we are interested 
in the conformal invariance inside the nontrivial region 2. In region 1 the critical point is 
gaussian, while it coincides with the SRI in region 3. The latter two cases are known to be 
conformally invariant, as will be reviewed in this and the next section. Understanding these 
arguments will be essential for the proof in the case of interest, in section 6. 


We will focus here on the proof of conformal invariance of the SRI critical point to all 
orders in the ^-expansion. Consider thus the usual, local theory in D = 4—e dimensions:^' 




d^x 




4! 



(4.1) 


Contrary to the initial conviction of some of us. 

^^We will use in this section not to confuse with 4i,d,e used in the sections dealing with the LRI. 

In the SRI case there is wavefunction renormalization, so we have to distinguish between the bare field yg 
and the renormalized y. 
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The model is regularized and renormalized using the dimensional regularization and the 
minimal subtraction scheme. There is an IR hxed point, called the Wilson-Fisher (WF) 
hxed point, which provides a held-theoretic description for the SRI critical point. We 
are specializing directly to the massless case, since in the considered scheme the mass 
renormalizes multiplicatively, and the hxed point is reached for zero bare mass. 

The idea of the proof [38, 39, 30, 40] is to construct the Ward identities for scale and 
conformal invariance broken by the effects of the running coupling, and to show that the 
breaking effects disappear at the hxed point. We will follow here a simplihed version of the 
presentation by Brown [30].^® 

Consider the canonical stress tensor of the theory: 

{\{dipof + (Ao/4!) ip^] . (4.2) 

It’s easy to compute its divergence and trace: 

= (4.3) 

T% = - (.0/2 - !)£ + (1/2 - D/4)9yS, (4.4) 

where the E and are operators proportional to the equations of motion of the theory: 

E = (fo + (Ao/3!)v3o} , E^, = d^(po {-^Vo + (Ao/3!)(/9^} . (4.5) 


Correlation functions of these operators are trivial: their insertions into iV-point functions 
of ip produce a bunch of ^-function at coincident points, see Brown’s (3.10), (3.28): 


E n 

6{x — Xi)G{xi... a; at) , 

i=l 

d 

G{xi ...xn; E^{x)) = 2^._^ <5(x - Xi)-^G(xi ... xjv) ■ 


(4.6) 


One may wonder why we bother at all about the operators E and E^, since they have 
vanishing correlation functions at non-coincident points. Indeed, in the usual CFT language, 
they would not even qualify to be called operators. However, here we are working in 
perturbative quantum held theory, and in this situation it turns out to be both legitimate 
and useful to have access to E and E^. Legitimate because in the regularized theory the 
equations (4.6) make perfect technical, and not just formal, sense. Useful because the fastest 
derivation of the Ward identities uses these operators, as we will see momentarily. 


We will need next Brown’s Eqs. (3.13), (3.14) and (3.18) which relate the nonrenormal- 
ized operator ipQ to the renormalized hnite operator 

^¥^0 = + (7(A)/£)E + const.d'^^pl . (4.7) 

^®Brown’s goal was to construct a finite, renormalized, stress tensor, which he then also used to discuss 
the Ward identities. Since here we are only interested in the Ward identities, we can derive them using 
the canonical stress tensor, differing from Brown’s finite tensor by total derivative terms. This allows us to 
avoid a complicated discussion of tensor operator renormalization, using only the rather simpler Eq. (4.7). 

^®In accord with the rest of this paper, we denote by /3(A) the full beta-function, equal to — eA-|-/3Brown(A). 
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Using this relation in (4.4), we get: 

T\ = ^;3(A)[/| - A{\)E + consf.dVi . (4.8) 

where A(A) = {D — 2)/2 + 7 (A) is the full (classical + anomalous) dimension of (p. 

One could (as Brown does) construct an improved stress tensor whose trace would not 
contain a term. This requires adding to the well-known improvement ieim.^{d^dy — 
At the classical level ^ = —(iA —2)/[4(iA —1)], but it gets corrected at the quantum 
level [30]. If one is only interested in the Ward identities, as we are here, all this extra work 
can be avoided, since will drop out. 

So, we proceed with the derivation of the Ward idenities. Consider the objects 

, C/ = T^u{2x'^x^ - S'^^x^) . (4.9) 

These are called scale and special conformal currents. By (4.3) their divergences can be 
expressed as: 


d^V^ = -x^E^ + T\, (4.10) 

= -{2x^^x^ - 5>^^x^)E^ + 2x^T^^^ . (4.11) 

Inserting these equation in an iV-point function and integrating over the full space, we get:^° 

J d^x G{xi ... a: at; -x^E^{x) + T^^{x)) = 0 , (4.12) 

f d^x G{xi... xn] —{2x^x^ — 5^^x^)E^{x) + 2x^T^^{x)) = 0 . (4.13) 


The contribution of in vanishes after integrating by parts. Contributions of the 
other terms can be all expressed using (4.8) and (4.6). We hnally obtain: 


A(A)]G(a;i.. .xat) =/?(A)^ f d^x G{xi... xn][p'^]{x)) , (4.14) 

i=i ■ 


E 



S‘‘P)-E+2A(\)xt 


G{xi. ..xn) = 2/3(A)^ 


d^x x^ G{xi.. .Xn', [p"^] (x)) . 

(4.15) 


These are the Ward identities expressing the breaking of scale and special conformal invari¬ 
ance by the running coupling. They are valid at all distances and for all A. For the purposes 

^°The boundary term arising when integrating by parts vanishes. To show this one has to estimate a 
correlation function of a group of ip’s and a widely separated T^^,. In the critical theory and for a traceless 
this would be l/|xp^ by using the stress tensor dimension. But here we need an estimate for a theory 
which has not yet reached the fixed point, and for an unimproved For a quick and dirty estimate, 

notice that there will be at least two propagators connecting to the y’s, with two derivatives acting on 
them. This gives decay ~ which is more than sufficient for D near 4. The part of gives 

an even smaller contribution. An analogous argument justifies integrating by parts a few lines below. 
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of studying the IR fixed point, we have to go to large distances and to take the limit /r —0, 
so that X{fi) —>■ A*. In this limit the beta-function multiplying the RHS vanishes. One is 
thus tempted to conclude that the IR hxed point is both scale and conformally invariance. 

While the conclusion is certainly correct, we believe that this last step of the argument 
merits somewhat more attention. There is a subtlety here. It’s true that /3{X) —)■ 0 as 
—)■ 0, but what about the integrals multiplying it? The integrals are //-dependent because 
the operator is normalized at the scale //. Could it be that the integrals grow in the 
/i —)■ 0 limit, overcoming the /3(A)/i^ suppression? This issue seems to have been neglected 
in the literature. In the next section we analyze the situation in detail and prove that this 
does not happen: the integrals behave as 0(/i“^), and so inferring the suppression from /9(A) 
is correct. Our argument is a bit technical, and it would be nice to hnd a simpler proof.^^ 

Taking that argument for granted, we have shown that the long-distance limit of the 
correlation functions of [p] transforms covariantly under the conformal group, with [p\ being 
a primary held of dimension A (A*). For a complete proof of conformal invariance of the 
SRI critical point, one needs to extend this result to correlators of composite operators. 
Correlation functions of composites can be related to those of the fundamental held by 
using the OPE, and it seems rather plausible that they will also be conformal. In appendix 
B we sketch a proof of this fact. 

At this point one may ask, analogously to what we did in section 3: can conformal 
invariance of the SRI critical point be tested by checking that two point functions of diherent 
primaries vanish? As mentioned in section 3.1, the correlator {p p'^) is nonzero in this theory 
because p'^ is a descendant of p. As for the {p"^ p^) correlator, we computed it and found 
that it does vanish at the IR hxed point to the hrst nontrivial order in e, which in this 
theory is This computation is a bit nontrivial as it requires careful analysis of the 

mixing between the operators p^ and d'^p'^, so we don’t present the details [43]. 

4.1 Scale and conformal invariance breaking effects 

Here we will study carefully how the broken Ward identities (4.14) and (4.15) imply scale 
and conformal symmetry of the WF hxed point. 


4.1.1 A proof which is too quick, and a paradox 

Let us hrst discuss a deceptively quick proof.Can’t we just set A = A* in these equations? 
On the one hand, the correlators in the LHS then become the IR limits of the correlators 

^^Such a proof can perhaps be given in the framework of the local CS equation [41, 42] for the theory 
with a;-dependent A coupling. We thank Hugh Osborn for sharing with us his thoughts and preliminary 
results in this direction. 

^^The argument given below equation (3.16) that 0(e) is trivial in the LRI does not apply in the SRI. 
The reason is that for e = 0 the SRI Feynman integrals are log-divergent. When moving to £ > 0 they have 
poles enhancing the effect. 

^^See e.g. p.29 of [40]. 
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of the full theory. On the other, the RHS side goes away as (3{K) = 0. Right? 

In our opinion, this argument is plagued with the following basic problem: the integrals 
multiplying /?(A) become UV divergent for A = A*. Indeed, for y ^ Xi we can use the OPE 

^{xi) X [ip%y) ~ I (4-16) 

\xi - y\ 

where A^4 = D + e + 0{e‘^) is the IR dimension of Thus the integrals diverge as A^, 
where A is the UV cutoff marking the boundary where the asymptotics (4.16) applies. In 

the strict A = A* limit we have A = oo, the RHS of the Ward identities becomes 0 x oo, and 

we learn nothing. 

Clearly, we should resolve this 0 x cx) ambiguity by keeping A slightly different from A*. 
We should then identify A with — the scale where the full RG flow transitions from the 
UV to the IR fixed point. This scale can be defined by the equation 

A(/ic)~|A,. (4.17) 

Since for \xi — y\ ^ the dimension of is near its UV fixed point dimension D — e, 
the UV divergence of the integrals will be cut off at Pc- We then send yc ^ oo, keeping y 
and X fixed. In this limit A approaches A* as 

K-\^e{y/ycY. (4.18) 

We obtain 

/d(A) X (integral) = 0(A* — A) x A^ = 0{y~^ x y^) = 0(1). (4.19) 

Here we kept track only of the yc dependence, as other parameters are kept hxed. 

We thus have a paradox: in the yc ^ oo limit the suppression from the /^-function seems 
to be exactly compensated by the growth of the integral. The RHS of the Ward identities 
seems to remain 0(1) and we show neither scale nor conformal invariance. 


4.1.2 Correct proof 

The resolution of the paradox is that the naive dimensional analysis does not capture the 
true asymptotics of the integrals multiplying the /3-function. The estimate 0{yl) applies 
separately to the UV and IR parts of the integral, but as we will now explain, a cancelation 
occurs when the two are added. The full integral is then 0(1) instead of 0{yl). 

From now on, instead of taking yc ^ oo with x, y fixed, it will be more convenient to 
keep yc fixed and approach the IR by taking x /:§> yc^, y yc- We will then have to keep 
track how various terms scale in the /i —)■ 0 limit. The basic issue—the insufficiency of the 
naive dimensional analysis—will however be the same as in the above paradox. 

Starting with the scale Ward identity (4.14), we notice that it can in fact be analyzed 
in two ways. One way is through the OPE as above: 

V?(0)b1(2/) ~ C{y, y)(p{0) . (4.20) 


20 



The argument /i reminds that [(^^] is normalized at the scale fi. In the a: —)■ oo limit the 
RHS of (4.14) can be approximated by using this OPE as 

7V/?(A)^ J d^yC{^^,y)xGix,...XN). (4.21) 

The second, actually more standard, way to analyze (4.14) is rewrite its RHS as^^ 

/3(A)(a/aA)G'(xi...a:^), (4.22) 

turning it into a CS equation, which can be solved similarly to the familiar two point function 
case. In particular this implies that at long distances the correlator is scale invariant with 
anomalous scaling dimensions: 

0(^2/!,..., ei/v) (X e » 1, (4.23) 

where the points yi,.. .yiy are assumed to have 0(/i“^) distances among them. 

Let us now take (4.23) and plug it into the LHS of (4.14). Assuming that all distances 
among x* are large compared to the LHS of (4.14) equals: 

Ar(A(A) - A(A*))G'(xi...x,v). (4.24) 

Comparing this to the RHS as given in (4.21), we get an interesting integral constraint on 
the OPE kernel: 

/3(A)L /'d“!/C(ft!;) = A(A)-A(A.). (4.25) 

Expanding this equation in A — A* for A(/i) near A* (i.e. p ^ pc), we get 

j d^yC{^i,y) = 0{1) (p < p^) • (4.26) 

It’s instructive to compare this result with simple dimensional analysis. By dimensional 
analysis, C(/i, y) must be of the form 

X f{\y\y, A). (4.27) 

Moreover, at large distances and to the leading order in epsilon it should scale oc l/\y\^^^, 
since the IR dimension of [(p^] is D + e + 0{e‘^). We will ignore the admixture of 
operator to [(p"^], since such admixture does not contribute to the integrals in the RHS of 
the Ward identities. We conclude that / oc (||/|/i)“^^ in the IR and thus 

C{ii,y) oc 1/[||/|^+V^], y > (4.28) 

Thus, were we to estimate the integral in (4.26) by extending the integration from to 
oo, we would get a result which is not 0(1) but grows in the /i —?■ 0 limit as: 

h" X = {yc/yY • (4.29) 

^^See [30], eq. (3.19); factor is missing in the RHS of that equation. 
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This is the same growth which plagued the argument in section 4.1.1. Apparently, there is 
a cancelation between the IR and the parts of (4.26). Separately they blow up in the 
/i —0 limit, but their sum remains 0(1). 

To summarize the logic so far: there are two ways to analyse the scale Ward identity. 
The hrst is based on using the OPE which reduces to the integral of the OPE kernel. One 
needs further input to get a good estimate for this integral. This input comes from using 
the second way of analysis, which is based on the CS equation and gives scale invariance. 
From consistency of the two ways we obtain (4.26). 

Now to the special conformal Ward identity. Its RHS cannot be reproduced by differen¬ 
tiating the Lagrangian.^® Thus there is only one way to proceed, that based on the OPE. 
We will now show, using (4.26), that the integral in the RHS of this Ward identity: 

^^yy^G{xi...XN][^%y)) (4.30) 

remains for /i —)■ 0 of the same order of magnitude as the quantities appearing in the LHS. 
Multiplied by /3(A) it vanishes, proving conformal invariance in the IR. 

We will choose all x* and their distances large, of the order s ^ We will also 

take /i /ic to have the coupling closed to its hxed point value. Finally, we will assume 
s /:§> /i“^. The rationale for this further separation of scales should become clear below. We 
then divide the integral into three regions: 

1. y is at distances < /i“^ from one ofxt. Here we can use the OPE. We have {z = y—Xi): 


if [ C{ii,z){z^ + x^) = Q + x^ X 0{l) (4.31) 

where we used that C{z) is spherically symmetric, and that the integral (4.26) remains 0(1) 
even when we truncate it to ?/ < Indeed, the truncated away part of the integral can 
be estimated by naive dimensional analysis using (4.28), and it’s 0(1).^^ 

2. y at distances between ii~^ and s. Here the OPE becomes less and less precise as the 
distance from Xi increases. But we can still use it as an order of magnitude estimate, by 
saying that the insertion of suppresses the correlation function by (see eq. (4.28)) 

0(l/[|.-r+y']). A.-'<k|<s. (4.32) 

This is enough to conclude that the integral over this region is of the same order of magnitude 
as the LHS. 

^®The integral converges in the UV since the UV dimension of f is D — e. 

^®Unless one introduces x-dependent couplings, see note 21. 

^^One could ask about the contributions of the subleading OPE terms, like those proportional to and 
to the other operators. All such terms can be estimated using the naive dimensional analysis and they are 
small. Indeed, they are power suppressed compared to the leading term by at least one power of y/s. Since 
the naive dimensional analysis was only missing the true behavior of the leading term by a power of e, it’s 
more than sufficient for the subleading terms. 
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3. y at distances 3> s from the group of Xi. The correlator decays as The 

integral is dominated by y near the lower integration limit, which is the same as region 2. 

This finishes the perturbative proof that the SRI correlations functions at the IR fixed 
point enjoy special conformal invariance. 

4.2 SRI conformal invariance beyond perturbation theory^^ 

Let us briefly discuss what is known about the conformal invariance of the SRI critical point 
beyond perturbation theory. First of all, in 2d the SRI is exactly solvable on the lattice. 
Polyakov [37] first noticed that the critical three point correlation functions following from 
the exact solution satisfy the restrictions of global conformal symmetry. In the continuum 
limit the 2d SRI critical point is a minimal model CFT [12]. Aspects of gobal conformal 
and Virasoro invariance of the SRI critical point observables were also proven rigorously by 
mathematicians, pioneered by S. Smirnov for nonlocal observables such as interfaces [44], 
and recently extended to correlation functions of local operators in bounded domains [45] . 

In 3d, the most dramatic evidence for the conformal invariance of the SRI critical point 
comes from its recent exploration using the conformal bootstrap methods [14-19, 13, 20- 
22]. These methods take conformal invariance as an assumption, which leads to a system 
of equations on the operator dimensions and OPE coefficients of the theory. The system is 
tightly constraining, and the fact that it has a solution which is in agreement with everything 
computed about the 3d Ising model using more pedestrian techniques is strong evidence for 
the validity of the assumption. 

Additional evidence for the conformal invariance of the 3d SRI critical point comes from 
the recent Monte-Carlo simulations of the model in the presence of a line defect [46] and in 
a ball geometry [47] . 

Let us now discuss prospects for a nonperturbative proof of the conformal invariance of 
the 3d SRI critical point. Because of the above evidence, we are certain that the model is 
conformally invariant. Still, it would be interesting to find a proof, because it may teach us 
about the conformal invariance of other models where the evidence is as yet lacking. 

Rather than deriving everything from scratch, it seems reasonable to allow a proof based 
on some assumptions. One rather minimal and natural set of assumptions is as follows:^® 

1. The 3d SRI critical point is a unitary quantum field theory, whose correlation functions 
satisfy the Wightman (or Osterwalder-Schrader in the Euclidean) axioms; 

2. It is scale invariant, with all local operators having well-defined scaling dimensions; 

3. It has a local stress tensor operator—a conserved rank 2 tensor operator of dimension 
exactly D = 3, whose correlation functions satisfy Ward identities. 

^®This section is outside of the main line of reasoning of this paper and can be left out on the first reading. 

^®In the future it would be interesting to prove the 3d conformal invariance without any assumptions. 
For the moment we follow the advice of Ennio De Giorgi: “Se non riesci a dimostrare I’enunciato, prendi 
una parte della tesi e mettila nelle ipotesi, e continua cosi finche non ci riesci.” 
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Because of assumption 2, the problem is often formulated as “Does scale invariance imply 
conformal invariance?”^*’ As lucidly explained long ago by Polchinski [52], this is not guar¬ 
anteed in this setup. The crucial question is whether the trace of the stress tensor contains a 
total derivative term in addition to terms vanishing at the critical point or at non-coincident 
points: 

T% ^ (?) (4-33) 

The operator is called a virial current. It must be a vector operator of scaling dimension 
exactly D — 1, to match the scaling dimension of The appearance of d^jV^ in does 
not contradict scale invariance, because the extra term vanishes when integrated over the 
full space, and does not disturb the scale Ward identity. However, the derivation of special 
conformal Ward identity requires integrating multiplied by see Eq. (4.13). The 
extra term then does not go away unless is itself a total derivative, 

From a modern perspective, then, to check if a theory is conformally invariant it’s 
enough to see if it contains a vector operator of dimension exactly D — 1 which is (a) a 
singlet under all global symmetries, (b) not conserved, and (c) not a total derivative. In 
weakly coupled examples it can be seen by inspection whether such an operator exists. As 
Polchinski [52] pointed out, the theory in D = 4 — e dimensions does not have any 
nontrivial candidate, and hence its IR fixed point must be conformally invariant to all 
orders in perturbation theory. 

In the strongly coupled situation one can argue (see e.g. [53-55]) that the (a,b,c) require¬ 
ments cannot generically be satisfied, and thus scale invariance generically implies conformal 
invariance. Indeed, what’s the likelihood that there will be a dimension D — 1 vector 
which is not conserved, given that all non-conserved vectors generically acquire anomalous 
dimensions? This being pretty obvious, the real question is to what extent one can get rid 
of the genericity assumption. 

Recently, the problem of showing conformal invariance of the 3d SRI critical point was 
revisited in [56]. Using an argument phrased in the language of the functional renormaliza¬ 
tion group, the authors arrive at the same conclusion as Polchinski—that scale invariance 
without conformal invariance is only possible if the theory contains a dimension D — 1 global 
singlet vector V. There is little doubt that their vector is the same as the virial current in 
Polchinski’s argument, although this connection is not made explicit in [56]. 

The rest of the discussion in [56] evolves around trying to show that a vector with 
dimension Ay = D — 1 is impossible. First they appeal to genericity. Compared to the above 
genericity argument, they add the following observation: to avoid conformal invariance, a 
dimension D — 1 vector should exist in the WF hxed point not only for D = 3, but in 
a whole interval of dimensions around 3, since otherwise one could argue conformality in 
D = 3 by continuity from D near 3. 

Second, they also point out that, independently of any theoretical arguments, the 
dimension of the lowest singlet vector can and should be measured in lattice Monte-Carlo 

^°See [48] for a review and in particular [49-51] for recent nontrivial progress in 4d. 

^^For example, the term d^ipQ in did not disturb the derivation of the conformal Ward identity in 
section 4 because it corresponds to Ri = d,i(yo)i ^ total derivative. 
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simulations. If it’s larger than D — 1, as it probably is, this will imply conformal invariance. 
Such a test of conformal invariance is probably easier, and more robust, than the checks of 
its predictions carried out in [46, 47]. 

To summarize, there is currently no nonperturbative proof of conformal invariance of 
the 3d SRI model critical point. Yet we are all but sure that it is conformally invariant. 

Note added: A third arXiv version of [56] just appeared which presents a proof of a 
lower bound Ay ^ D — 1 + t], based on Lebowitz inequalities. If conhrmed true, this would 
constitute a rigorous nonperturbative proof of conformal invariance in the 3d SRI. 


5 Conformal invariance of the gaussian phase 

The nonlocal gaussian theory described by the action (1.2) is also conformal. This is widely 
known in high energy physics,and more recently has also been discussed from statistical 
physics perspective. In this section we will review this fact pedagogically, trying to bridge 
the gap between the two communities. 


5.1 Direct argument 

We will start in the antichronological order. As pointed out in [57] , the fractional Laplacian, 
just like the ordinary Laplacian, is covariant under conformal transformations. Namely if 
X —)■ x' is a conformal transformation and 

(j)'{x') = \dx'/dx\~^'^^'^(j){x) (5.1) 

with Afp given in (1.3), then 

C!^^[x!) = \dx'La^[x). (5.2) 

As a result the action (1.2) is invariant under conformal transformations. The proof of (5.2) 
is given in [57] and we will not repeat it here. As usual, covariance under translations, 
rotations, and dilatations is obvious, and a simple calculation establishes covariance under 
the inversion. 

Notice that since we are dealing with a nonlocal theory, there is no reason to expect 
that in d = 2 the global conformal invariance of the action (1.2) gets enhanced to the full 
Virasoro invariance, and indeed this does not happen. 

^^There, this theory is sometimes referred to as Mean Field Theory or Generalized Free Field. 

Although it’s not essential for this paper, we would like to point out that there are also examples of local 
2d theories which have global conformal but not Virasoro invariance. Once such theory is the “biharmonic 
scalar” with the Lagrangian Its stress tensor trace is of the form which is enough for global 

conformal invariance but not enough for an improvement to make it traceless and recover full Virasoro, since 
in 2d one needs = d^Y for the latter [52]. See [58, 57] for a discussion. 
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5.2 Argument from correlation functions 

The previous argument shows the invariance of the action. Since the theory is gaussian, 
conformal covariance of the correlation functions follows. It’s also easy to check the trans¬ 
formation properties of the correlation functions directly. This way of proving conformal 
invariance predates the previous one. One of us (S.R.) was aware of it at the time of writing 
[59], and an analogous discussion for a nonlocal vector theory appeared in [60]. 

Let’s start with the two point function of 0. It is given by |a; — which indeed has 

the form of a two point function of a primary scalar of dimension A^. The iV-point functions 
of (j) are given by Wick’s theorem, since the theory is gaussian. A moment’s thought shows 
that since the two point function transforms as it should, the A-point functions will do so 
as well. So all correlation functions of 0 are consistent with conformal symmetry. 

The theory contains more operators, for example the normal ordered products 
Their correlation functions are dehned by just leaving out the Wick contractions at coinci¬ 
dent points, hence they will also be conformally covariant. 

Although there are still more operators in addition to the ones considered above, all of 
them can be obtained by taking repeatedly the OPE of 0 with itself. Correlation functions 
of these operators will inherit conformal transformation properties from the correlators of 
0.^^ Hence the theory is conformal. 

5.3 AdS argument 

Going further back in time [61, 62], an easy way to see conformal invariance of the theory 
(1.2) is from its dual holographic description. This description is given by the free massive 
scalar held T = 

^Ads = ljd'^x + M2VI/2] (5.3) 

in the d + 1 dimensional AdS space with the additional coordinate z > 0 and the metric^^ 

ds^ = z~‘^{dx^dx^ + dz^). (5.4) 

The dual theory (5.3) is manifestly invariant under the conformal group SO{d -|- 1,1). This 
is because AdS is a symmetric space, and SO{d + 1,1) is its group of isometries. Since the 
action (5.3) is written using the generally covariant rules, isometries of the underlying space 
become symmetries of the action. 

To show that the theory (5.3) is indeed dual, i.e. equivalent, to the nonlocal theory (1.2), 
we must “integrate out” the bulk of the AdS space by solving the equations of motion of 

^^See the argument in appendix B for the WF fixed point. 

^^This metric is considered fixed, rather than fluctuating like in the conventional AdS/CFT correspon¬ 
dence. This is related to the fact that the gaussian theory does not have a local stress tensor operator. 
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the field z), and derive an effective action for its bonndary valne 0b at 2 : = 0. The held 
T(a;, z) behaves near the bonndary as 

T(a;, z) ~ z‘^~^(j)b{x ), (5.5) 

where A is a root of the qnadratic eqnation A(A — d) = M^, 

A± = l{d±Vd^ + 4M^). (5.6) 

Picking A = A+, the effective action has the form identical to (1.2) [61, 62] 

S'eff oc J d'^xd%(j)b{x)(j)b{y)/\x. (5.7) 

The range of dimensions A is however not qnite the same. In fact should satisfy the so- 
called Breitenlohner-Freedman bound > —d‘^/4 [63], which allows to reproduce all the 
dimensions above d/2, while in (1.2) we need dimensions in the range (d — 2)/2 < A < d/2. 
As explained in [64], dimensions in this range are reproduced by taking the smallest root 
A = A_, which in this setup requires adding a boundary term in the action (5.3). 

Since the original work of [61, 62], the above way of thinking about why the nonlocal 
gaussian scalar is conformal has been used numerous times in the high energy physics 
literature. See e.g. [59, 65-67] for a partial list of recent references. This theory is very well 
understood from the CFT perspective, e.g. in [65, 67] the conformal block decomposition 
and the 0x0 OPE have been worked out explicitly for any A. 

5.4 Caffarelli-Silvestre trick 

The idea of the previous argument is to rewrite the nonlocal gaussian theory as an equivalent 
higher dimensional theory whose conformal invariance is manifest. We will now present a 
second way to implement this idea, based on an observation of Caffarelli and Silvestre [68]. 

Consider a scalar field $ = <F(a:, y) where the extra coordinate y now takes values in the 
fiat Euclidean space of p = 2 — a dimensions. Hopefully the reader is not disturbed by the 
fact that p is in general fractional. In this space we consider the massless scalar field action: 

Acs = j d^x d^y [(0,<h)2 + [dy^f] . (5.8) 

We will now show that this local action is equivalent to the nonlocal action (1.2). 

Let (/{x) be the value of the field <F on the d dimensional hyperplane p = 0: 

0(a;) = <F(a;, 0), (5.9) 

and consider the effective action for 0 obtained by integrating out the rest of the space. 

^®This is the nonlocal part of the effective action. As is well known, the on-shell effective action contains 
divergent local counterterms which are implicitly subtracted here. 
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To find it we have to solve the equations of motion of the higher-dimensional theory 
subject to the boundary condition (5.9). It’s clear that we can restrict to the sector of helds 
radially symmetric in the y variable. In this sector the action becomes [z = \y\ > 0): 

Scs = Spj + (5.10) 

Recall that Sp is the unit sphere volume in p dimensions. The equation of motion is 

dl^ + {l-a)d,^/z + dl^ = {). (5.11) 

The solution decreasing at large 2 ; takes in momentum space the form; 

<I)(p, z) = const.{\p\zY/‘^K„/ 2 {\p\z) 0(p). (5.12) 

For small z this has an expansion: 

<h(p, z) = (1 + const.{\p\zY -\-.. .)4>{p). (5.13) 

The ... terms are 0{zY compared to the shown ones and are subdominant in the range of 
interest 0 < cr < 2. We integrate by parts in (5.8), pick up the boundary term, and hnd 

Scs = -^pj j d^p4>{p)\p\‘"4>{-p), (5.14) 

which is the nonlocal action (1.2) in momentum space. 

It is now easy to complete the argument for conformal invariance. We started with a 
massless scalar theory in flat space (5.8), which is conformally invariant in an arbitrary 
number of dimensions. When we construct an effective theory by integrating out the space 
away from the y = 0 hyperplane, we are guaranteed to obtain a theory invariant under the 
subgroup of the d + p dimensional conformal group which leaves invariant this hyperplane. 
This is precisely the d dimensional conformal group. 

It’s interesting to compare this argument with the AdS one in the previous section. The 
effective action computation here and in the AdS case proceed rather similarly. In fact 
the equations of motion of the AdS theory (5.3) can be mapped onto those of (5.10) by 
dehning <I)(a:, z) = z~^’i=d/{x, z). A minor difference is that the on-shell action is hnite in the 
Caffarelli-Silvestre case, while in the AdS case it needs to be renormalized with boundary 
counterterms which we implicitly subtracted in (5.7). 

There is however a big difference in explaining why the equivalent higher-dimensional 
theory is conformally invariant. In the AdS case we appealed to the isometries of the 
underlying geometry, while in the CS description we use the fact that the massless scalar 
theory in flat space is conformally invariant. 


^^The form of this effective action could also be predicted from the following argument. The free scalar 
<i> two point function depends on the distance as oc The same dependence is inherited by (p when 

we set y = 0, and we recover the expected two point function of the nonlocal scalar field of dimension d — a. 



Concerning the prior history of the given argument, Caffarelli and Silvestre noticed 
that the equation of motion (5.11) gives rise to the fractional Laplacian and showed the 
equivalence of the two actions. They also noticed that (5.11) is nothing but the Laplace 
equations for functions radially symmetric in p extra coordinates. The Caffarelli-Silvestre 
description in the radially reduced sector was used by Rajabpour [57] to discuss some aspects 
of the nonlocal gaussian theory. Using the full d + p dimensional description to argue for 
the conformal invariance of the theory seems to be done here for the hrst time. 


6 Conformal invariance of the LRI critical point 

We will now present a proof that the LRI critical point is conformally invariant also in the 
nontrivial region 2. After all the preparatory work in sections 4 and 5, the proof looks 
almost inevitable. 

The key idea is to rewrite our theory (2.3) as a defect quantum field theory. Namely, 
according to the discussion in section 5.4, we can rewrite the action as^® 

S=lJ + ^ ( 6 . 1 ) 

In the hrst term <l> is the free scalar Caffarelli-Silvestre held from section 5.4, dehned on the 
d = d + p dimensional space X = [x, y). The second term represents interaction living on 
the defect: the d-dimensional plane located at y = 0. The number of extra dimensions will 
be fractional in the case of interest (a near d/2), but in spite of this fact the theory (6.1) 
makes perfect sense, at least in perturbation theory.^® 

In the theory (6.1) we will consider V-point functions of *h(W), G{Xi... X^). The 
correlators of the original theory (2.3) can be obtained by taking the ?/ —t 0 limit. 

The hrst part of the proof will be to derive broken scale and conformal Ward identities 
satished by the correlators G{Xi.. .Xj^). We will then discuss how these Ward identities 
imply conformal invariance of the IR hxed point. 

6.1 Ward identities 

Crucially, the theory (6.1) is local and so its Ward identities can be derived by a variation 
of the method from section 4. We start by considering the canonical stress tensor: 

Tmn = . (6.2) 

The indices M, N... will run over the full d-dimensional space, p,i>... over the d-dimensional 
“parallel” subspace, and m,n... over the p-dimensional “perpendicular” subspace. The 

^^Normalization of the first term is different from (2.3) where it was fixed via (2.4). This difference is 
unimportant for the proof of conformal invariance. 

^®See section 6.4 for a nonperturbative discussion. An introduction to perturbative quantum field theory 
in non-integer dimensions can be found in [24] . 
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is the Kronecker delta in the parallel subspace: Sj^j^ = if both indices are parellel, and 
zero otherwise. 

Unlike in section 4, here we don’t make a distinction between the bare and renormalized 
held <h, because does not get renormalized: since the interaction term is located on the 
defect, it cannot renormalize the bulk kinetic term. The boundary kinetic term corresponds 
to an irrelevant operator and cannot be generated either. 

The next step is to hnd the divergence and the trace of Tmn- Direct computation gives: 

O^^Tmn = -En + , (6.3) 

T^m = + (1/2 - d/A)dl ^^, (6.4) 

where we recognize the operators proportional to the equations of motion: 

E = $ {-a|4) + ^ ^ + 5(^’)(|/)(^/o/3!)4)3} , (6.5) 

and a new object called the displacement operator. It is given by: 

= (W3!)<l>'9n<l> if A = n (6.6) 


is a perpendicular index, and vanishes otherwise. This operator represents an inhnitesimal 
movement of the defect in an orthogonal direction. 

Unlike <h, the operator *h'^|y=o = 0“^ appearing in (6.4) is not a hnite operator; it will be 
related to a hnite renormalized operator [ 0 "^] via a rescaling 

0" = ^4[0l . (6.7) 


Contrary to (4.7), in the LRI there is no other operator with which 0^ could mix. The 
coefficient Z 4 is easy to compute. For any hnite correlation function we have 


A 

dg 




(...) 

=const 

—(•••) 

dgo^ ^ 


(■ 

(• 


d‘^a:/i"[ 0 ^](a;)), 
d'^a; 0 ^(a;)). 


( 6 . 8 ) 


On the other hand, diherentiating (2.5) with respect to log /i we obtain the relation dg^/dg = 
—ego/(3{g). It follows that: 

Z 4 =-0(^)p7(e^o), (6.9) 

which is identical to the hrst coefficient in (4.7). Plugging this into (6.4) we obtain: 

T^m = -A^E + (0(^)/4!)/i^7^)(7/)[0"] + (1/2 - J/4)0|<F2 . (6.10) 


We next consider the dilatation and special conformal currents: 


rji \rN /7 L rji f(r\\rN \rL ^NL\r2\ 

— J-MN^ 5 ^ —0 A J. 


( 6 . 11 ) 
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The divergence of the scale current is given by 

d^VM = -X^{En + S^P\y)DM) + T^m = -X^Em + m • ( 6 . 12 ) 

The term proportional to Hat is seen to vanish, since either is a parallel index and then 
Dn = 0, or else it’s a perpendicular index and then X^5^\y) = 0. The divergence equation 
takes the same form as (4.10), and we obtain the scale Ward identity in full analogy to (4.14): 

N ^ „ 

Y,[X^.^x, + A^]G{X,...XN)=/3{g)^ / G{X,... Xx,[<p%x)). (6.13) 

i=i ■ 

Analogously let’s analyze the divergence of the special conformal current 

d^CM^ = -{2X^X^ - 6 ^^X^){En + 6^\y)DN) + 2 X^T^m ■ (6.14) 

We would like to derive a Ward identity corresponding to special conformal transformations 
leaving the defect invariant, i.e. when L = A is a parallel index. The extra term proportional 
to Dn then drops out just as for the scale current, since Dx is nontrivial only if iV = n is 
a perpendicular index. Then 5^^ = 0, while X^X^ vanishes when multiplied by 5^^\y). We 
therefore have the same divergence equation as (4.11), and an analogous Ward identity:^'^ 


E 

i=l 

= Jd^xx^GiX,...Xx;[(l>^]{x)). (6.15) 


(2X 


Mxx _ 5 


M\ 




G(A,...A 


N 


Now, the correlators G{Xi... Xn) behave continuously in the limit ?/ —)■ 0. This can be 
seen diagram by diagram in perturbation theory. It’s also closely related to the fact that 
<h does not acquire anomalous dimension, and thus its bulk-to-defect OPE is non-singular 
(see a related discussion in section 6.3 below). 

Thus, if one is primarily interested in the correlators at the defect, which are the 
correlators of the original theory, then one may set j/ = 0 in (6.13) and (6.15), obtaining 
restricted Ward identities satished by the correlators of the original theory (2.3): 


Y^[xi.d^.+A^]G{xi.. .Xn) = I d'^xGixi.. .xn-,[(P%x)) , (6.16) 

N 

E 

(6.17) 


d 


{2xfx^ - + 2A^x] 




G{xi.. .Xn) = 213(g)— / d'^x x^ G{xi... XN-,[(j)]{x)) 


At this point, if one wishes, one may altogether forget about the construct of the extra 
dimensions. Notice however that without this construct, it would remain pretty mysterious 

^°It’s not hard to convince oneself, like we did in section 4, that various boundary terms appearing when 
integrating by parts in the derivation of the Ward identities vanish by a good margin. 
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why the nonlocal LRI theory should satisfy Ward identities which are almost identical in 
form to the ones satished by the local SRld^ 

6.2 Prom the Ward identities to conformal invariance 

We will now discuss how the broken Ward identities imply that the IR hxed point is scale 
and conformally invariant. Naively this follows from the fact that P{g) multiplying the 
breaking terms vanishes at the hxed point. However, as mentioned in section 4 for the SRI 
case, this is too quick. A careful argument must show that the integrals multiplying (]{g) 
do not overcome the suppression. For the SRI, this was done in detail in section 4.1, and a 
modihcation of that argument will work for the LRI. 

Let us hrst consider the y = 0 correlators. The required modihcations are then minor. 
The proof proceeds via the OPE 0x [0“^] D 0. Using the scale Ward identity (6.16), rewritten 
as a CS equation, one derives scale invariance, and also obtains an integral constraint on 
the OPE kernel. There is a simplihcation due to the fact that does not depend on g 
in the LRI. Thus the LRI analogue of Eq. (4.26) will have 0 rather than 0(1) in the RHS. 
One then uses the integral constraint to show that the RHS of the special conformal Ward 
identity (6.17) is small in the IR. This part of the argument is identical to that for the SRI. 

As for the correlators G{Xi... Xjy) with helds inserted away from the defect, it seems 
a bit awkward to generalize the SRI argument to this case.^^ Fortunately there is an easier 
alternative argument: conformal invariance of these correlators can be seen to follow from 
two facts: (1) that for y > 0 they satisfy the (conformally invariant) Laplace equation, and 
(2) that their boundary values ai y = 0 are conformally invariant, as we already proved. 

Furthermore, conformal invariance of the composite operator correlators should now 
follow from by considerations based on the OPE. This is especially obvious for the bulk 
correlators, since one then needs only the bulk-to-bulk OPE, identical to that of the free 
theory. Composite operator correlators on the defect can be obtained from those away from 
the defect via the bulk-to defect OPE (see the next section). Alternatively, one can access 
such correlators via the OPE of fundamental helds on the defect. This last way of arguing 
would be completely analogous to the one for the SRI given in appendix B. 


6.3 Bulk-to-defect OPE and the anomalous dimensions 

Up to now we were using the theory (6.1) as a formal vehicle for proving things about the 
LRI. We would like to demonstrate, by a simple example, that this formulation can also be 

is possible that the Ward identities (6.16), (6.17) can be proved by endowing the LRI theory with 
some sort of nonlocal stress tensor operator. See [57] for steps in this direction for the gaussian case. Still, 
the most natural path to the nonlocal stress tensor lies through the Caffarelli-Silvestre construction. 

^^One would have to consider an OPE of the operator on the defect with ^{x,y) in the bulk. The 
cases of y small or large compared to and the distances among Xi appear to require separate treatments. 

^^This argument should be applied to the connected part of the bulk correlators, which satisfy the Laplace 
equation in each variable X^ independently, without contact terms in the bulk. 
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used for concrete computations. 


One additional concept which will be necessary in this regard is the bulk-to-defect OPE.'^'^ 
Let C>(x, y) be a dimension A bulk operator of the theory. The bulk-to-defect OPE says that 
inside correlation functions insertions of this operator can be replaced by an inhnite series 
of operators Ok inserted at the defect point x. By dimensional analysis, this expansion will 
take the form 

= (6.18) 

where the index (0) marks bare defect operators and their classical dimensions. For simplic¬ 
ity we only consider scalar operators and we will ignore operator mixing. The renormalized 
bulk-to-defect OPE takes the form 

k bl 

with the renormalized defect operators dehned as 

= Zk[Ok], ( 6 . 20 ) 

whereas the bulk operators do not renormalize. Clearly, then 

9) = 'ik{g)f^{yy, g ), (6.21) 

which ties the anomalous dimension of the defect operators to the bulk-to-defect OPE. At 
the IR hxed point the bulk-to-defect OPE will take the form: 

y) = Y^ [Ok] (x), (6.22) 

k \y\ 

This OPE will govern the asymptotics of any correlation function of 0{x,y) in the limit 
when y is much smaller than the x distance between the operators. 

As an example we consider the two point function of with both operators inserted 
away from the defect. At 0{g) this correlation function is given by 

+ . (6.23) 

= + + (6.24) 

The I integral is hnite for yi ^2 > 0. We can evaluate it using Feynman parameters: 

/[xi 2 , 1 / 1 , 2 / 2 ] = - —[m(1 - u)x^ + uyl + (1 - u)yl] 

(6.25) 

^^For a similar concept of bulk-to-boundary OPE in boundary CFTs, see for example [69, 70]. 
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To simplify, let’s consider the special case yi = 1/2 = y- We Taylor expand in perform 
the u integral for each term, and then resnm the terms in the Taylor series. This leads to 

r(4Al)~ , (6,26) 

We are interested in the limit y 0 : 12 . Becanse everything is hnite we may set e = 0, so 
= (i/4. For small y we hnd then that 


^[^i 2 ,y,y] = -Srf^-b + ^(c^/2) + log(|/Va^? 2 )] + 0 ( 2 /^). (6.27) 

^12 


We can snbstitnte this resnlt in the fnll correlation fnnction and compare this with the 
behavior expected from the bnlk-to-defect OPE: 


(^^(xi,y)^^(x2,y)} 


—)■ 


a 


y2A-2A 


X 


2A ■ 
12 


(6.28) 


Here A = nA^, while A is the dimension of the hrst defect operator, i.e. [(/>"■]. To the leading 
order in e that we are working, (6.27) shonld match (6.28) expanded in the limit of small 
anomalons dimension 7 „ = A — A. We also have to set g ^ in (6.23). We obtain: 


Cn — n! [1 + g^:Kn{'y + 'ip{d/2)) + 0(e^)] , 7 n — g*Kn + O(e^), (6.29) 


where Kn is the same as in (2.17). The anomalons dimension resnlt agrees with (2.18). It 
is cnrions to see the nsnally scheme-dependent terms 7 and 'i/((i/ 2 ) appear in the sqnared 
bnlk-to-defect OPE coefficient Cn- 


6.4 Beyond perturbation theory 

Finally, we wonld like to discnss, very briefly and in parallel with section 4.2 for the SRI, 
the prospects of a nonpertnrbative proof of conformal invariance of the LRI critical point. 

First of all one wonld have to be convinced that the extra-dimensional defect QFT 
formnlation ( 6 . 1 ) makes sense beyond pertnrbation theory. 

Alternatively, one can try to get rid of the fractional dimensions altogether by restricting 
the 2 / > 0 theory to the radially symmetric sector described by the rednced action (5.10), 
which is (i -|- 1 dimensional. This action shonld still contain the conformal symmetry of 
the original action. The disadvantage of this formulation is the explicit dependence of the 
action on y. On the other hand it’s manifestly nonperturbatively well dehned. 

Whatever the formulation one uses, one will be reduced to studying the Ward identities 
for the corresponding stress tensor operator. Analogously to section 4.2, one can ask: what 
if the stress tensor trace contains the following term: 

M ^ 5^^\y)d^v>^ (?), (6.30) 

^®For a collection of recent work trying to make nonpertnrbative sense of various aspects of quantum field 
theory in fractional number of dimensions see [2, 71, 72]. 
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where is a dimension d — 1 vector operator on the defect—a defect virial current. By 
assumption, is not identically conserved. Like in the discussion following the analogous 
Eq. (4.33), we see that unless is a total derivative its presence precludes conformal 
invariance of the correlation functions, while preserving their scale invariance. 

So, can our theory contain a defect vector which has dimension exactly d — 1 and is not 
a total derivative? Generically, this appears unlikely (with or without the total derivative 
clause). In perturbation theory, there was no candidate for such an operator (and so it’s 
not surprising that no such term is visible in (6.10)). To establish this rigorously and 
nonperturbatively appears as hard as the corresponding problem for the SRI. 

We hnally note that in the context of boundary rather than defect quantum held theory 
the question whether scale invariance implies conformality was recently discussed in [73]. 


7 Conclusions 

The general problem of scale versus conformal invariance is continuing to provide food for 
thought. In this paper we have discussed this issue for the specihc theory corresponding to 
the LRI at criticality. Since this model most notably does not have a local stress tensor, 
even the standard (genericity) arguments are invalid, leaving the question of scale versus 
conformal invariance hanging in midair. 

In the hrst part of the paper we have provided nontrivial evidence for conformal in¬ 
variance of the LRI by showing that the and (0^0^) two point functions vanish at 

criticality, at least up to terms of 0{e^). In the second part we were able to prove conformal 
invariance at the critical point to all orders in perturbation theory. The salient point of 
the proof was the construction of the LRI model as a defect theory in an auxiliary higher- 
dimensional space; this allowed us to work with a higher-dimensional stress tensor and 
construct a proof analogous to the one used for the Wilson-Fisher hxed point. Much like 
the SRI, it is plausible that the LRI is also nonperturbatively conformally invariant. 

What does our proof teach us about general long-range lattice models? First of all, 
the models that have an epsilon expansion starting from a generalized free theory, like the 
long-range 0{N) or Potts models,^® will be conformally invariant at their critical points to 
all orders by simple extensions of our proof. 

In hindsight it is now also easy to cook up long-range models which will not be confor¬ 
mally invariant. The simplest one is based on a vector rather than a scalar. The gaussian 
long-range action then has two terms consistent with rotation invariance. Equivalently, the 

^®The q-component Potts model has a Landau-Ginzburg description via a {q — l)-component scalar field 
with quartic couplings preserving Sq symmetry [74]. This description is trivially generalized to the long- 
range case. The phase transition is second order if the 5'q-invariant cubic operator is irrelevant at the IR 
fixed point. 

■^^Recall once again that in 1 -|-1 dimensions conformal invariance without a stress tensor implies a global 
SO{3, 1) invariance, not the full Virasoro symmetry. 
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vector two point function also has two terms: 

= {a5^y + hx^x''/x^)/\x\^^’^. (7.1) 

This agrees with the two point function of a conformal primary vector for h/a = —2, and 
with that of a derivative of a scalar primary for h/a = —2 Ab- For all the other ratios this 
gaussian model will not be conformal [60].'^® Let us now perturb the gaussian theory by a 
quartic local term {B^B^Y flow to the IR hxed point. We conjecture that this fixed 
point will be an interacting scale invariant theory without conformal invariance. It would be 
interesting to check this by explicit epsilon expansion computations of two point functions, 
like we did in section 3. 

In recent years the use of the conformal bootstrap has emerged as an excellent tool to 
analyze CFTs in various dimensions [14]. The main ingredients in this approach are the 
conformal block decomposition of four point functions and the associated crossing symmetry 
equations, which all follow directly from conformal invariance (see e.g. [54]). The present 
work therefore opens the door towards an analysis of the LRI using these methods. Such 
an analysis would be completely non-perturbative, and is likely to result in a high-precision 
determination of its critical exponents. 

With an eye towards the conformal bootstrap let us finally discuss the critical point of 
the LRI from the perspective of an ordinary CFT. Firstly, the LRI distinguishes itself from 
the more familiar, local, CFTs by the absence of a local stress tensor operator. Secondly, 
we expect a one-dimensional family of solutions to the crossing symmetry, parametrized by 
the exactly known scaling dimension of 0. Thirdly, our analysis uncovered an interesting 
fact that the second Z 2 odd primary operator is related to 0 through the nonlocal EOM: 
0^ oc €„</). This hxes the scaling dimension of to be + a. Curiously, this can be 
equivalently rewritten as 

A^3 = d-A^. (7.2) 

This identihes as the so-called ‘shadow’ operator of 0, which by dehnition transforms in 
a conformal representation with the same value of the quadratic Casimir as 0. The shadow 
operators often appear in discussions of conformal blocks (see e.g. [76]), but usually as a 
formal tool, since they do not belong to the local operators of the theory. On the contrary, 
Eq. (7.2) means that both 0 and its shadow are good local operators of the LRI critical 
point. 

One may wonder if the nonlocal EOM also implies an exact proportionality relation 
between the coefficients with which 0 and 0^ appear in any OPE. Unfortunately, this appears 
not to be the case. Using the nonlocal EOM once, we can easily relate the three point 

■^^One can also reproduce this theory a la Caffarelli-Silvestre from an auxiliary vector theory in a higher¬ 
dimensional space. The auxiliary theory action has two terms and (OmBm)^■ As is well known, 

such a vector theory without gauge invariance (“theory of elasticity”) is not conformally invariant [75, 53]. 

^®Other examples of theories with pairs of operators satisfying the “shadow relation” Ai -|- A 2 = d can 
be found among SUSY CFTs. Namely, such a relation can emerge if Oi is a chiral primary of dimension 
d/n where n is an integer, and O 2 = ((bi)"“^. We thank Leonardo Rastelli for sending us a list of such 
theories. One simple example is the W = 2 3d Ising model where n = 3. What sets LRI apart is that the 
shadow relation appears as a consequence of the nonlocal EOM. 
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functions of 0^ to those of 0. However, we then have to relate the normalizations of the 
0 and 0^ two point functions. This requires using the nonlocal EOM twice, subtracting 
the unperturbed 0 two point function as in (3.37). The hnal answer then depends on the 
relative normalization of (00) and (00)gg, which is unknown in general. See appendix C. 

We conclude with a brief comment about the limit a —?• a*. What happens with the 
nonlocal CFT describing the LRI critical point when we approach this limit from below? 
As we discussed in the introductions, for a > a* the LRI critical point is supposed to be 
in the SRI universality class. It would be simplest if all correlation functions continuously 
transitioned to the SRI ones in the limit. However, this seems problematic in view of the 
presence of 0^ as a primary in the LRI CFT. 

Consider for example the correlation function (0 0^ 0 0^) . In the LRI, its conformal block 
decomposition is expected to contain the contributions of two relevant Z 2 odd scalars: 0 
and 0^. On the other hand, is a descendant in the CFT describing the WF hxed point, 
while the second Z 2 odd primary operator (ip^) is irrelevant. It could be that 0^ decouples 
in the cr —?■ a* limit, but based on the discussion in appendix C this also seems unlikely. We 
thus are led to conclude that there is a discontinuity in the transition from the LRI to the 
SRI correlation functions. It would be very interesting to investigate this in more detail. 
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A The hard diagram 


In this appendix we will discuss the computation of the coefficient P 2 C in (3.43), which comes 
from the last diagram in (3.41). We need to evaluate the integral 



1 

1 

1 

1 

1 

T 

\x — 

\x — z\ 

p" \y — z\ 

|2a |^|2/3 

0|2/3 

2; 3d/2-7e/2 


where a = = {d — e)/4, 0 = 2 q;. We will show that 

T = 47r‘‘r(-^)/r(f)+0(€). (A.2) 

To get from here to P 2 C one needs to multiply by the symmetry factor 2^3^/(4!)^ = 6. 
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We will present two methods to compute T: one using Mellin-Barnes representations 
and another expansions into Gegenbauer polynomials. 


A.l Mellin-Barnes method 


The y integral in (A.l) has the following Mellin-Barnes (MB) representation ([77], Eq. (24)): 


where 




1 


1 


X 


d 

TTZ 


2 /P" \y-z 
dsdt 


2a 


\y\2y 


\Z — X 


\d-2e 


M{s, t) 


(27rz)2 r(e)r(^) r(^) 



(A.3) 


Ai(s,t) = r(-s)r(-t)r{-^ + ^-s 


^ , d 3e 


r I s -|-1 -|- 


d — e 


r I s -|-1 -|- 


d-2e 


(A.4) 


One way to get this representation is to introduce an extra factor \y — into the 

integral ly. The resulting integral ly is then conformal as the exponents sum to 2d. Its 
MB representation is given by Symanzik’s formula (Eq. (5) in [78]). Eq. (A.3) follows by 
taking the limit tc —)■ oo in which the conformal cross-ratios in Symanzik’s formula tend to 
|a;p /|;2 — xp and \z\^/\z — xp. 


The contours in (A.3) run from —ioo to ioo separating poles of the T functions whose 
arguments involve positive and negative s, t. We can take straight-line contours s = Sq + 
and t = to + iM., where (so,to) lies in the small red triangle in Fig. 2. E.g. Sq = to = 
(“f + t) ~ I ^dl do. That straight line contours can be chosen is related to the fact that 
the integral ly is both UV and IR convergent and does not require analytic continuation. 


We must now compute the integral over z, which is of the form 


/d-^- 



1 

ttI r(Ai + A2-|)r(f-Ai)r(|-A2) 


2Ai 

z — x\ 

p^2 r(Ai)r(A 2 ) 

\x\ 

| 2 (Ai+A 2 -f) r(2h-Ai-A2) 


with 

Ai = (d — e)/2 — f , A 2 = 3d/4 — 5e/4-|-s-f f . (A.6) 

However, there is an issue. In order to commute the integral over z with the MB integral, 
we must make sure that the ^ integral converges in the ordinary sense for all values of s, t 
in the integration contour. The integral converges for 


0 < Re(Ai),Re(A 2 ) < d/2, Re(Ai + A 2 ) > d/2 , (A.7) 

while for the other parameters it’s dehned by analytic continuation. This condition is 
violated by the above choice of contour. The solution is to shift the contour before doing 
the z integral. One possible shift is indicated in Fig. 2. The t contour is shifted up so that 
it crosses the horizontal line of poles at t = ^ and f = 0 and ends up in the bulk of 

the allowed grey triangular region. The new to satishes the condition 0 < fo < -f +I Ai¬ 


ds 
















Figure 2: The horizontal, vertical, and diagonal straight lines show the location of the poles 
in (A.4). In this plot d = 3 and e = 0.3, but the same picture is valid for all d of interest and 
e <C 1. Only the leading poles are shown; each line is accompanied by a series of parallel lines 
spaced by 1 in the direction shown by short arrows. The real parts (so)^o) of the integration 
contours in the MB representation (A.3) must lie in the small red triangle. On the other hand 
{so,to) must lie in the large grey triangle for (A.7) to be satisfied. The dashed arrow shows a 
possible contour shift, crossing two horizontal pole lines at t = — | + ^ and t = 0. 


After this contour shift, we perform the z integration and obtain 

T =--—(/i + I2 + h), 


r{e)r i^-f Y T i^-f ) 


where 


h = 
h = - 
h = - 




dsdt 
(27ri)^ 

-Res.__d , t ), 

27rz 2 

^Rest=oM(s,t) , 


M(s, t) = A4(s, t) X 


r(f-; + >.)r(| + f)r(-| + |-.-f) 


(A.8) 

(A.9) 
(A. 10) 
(A. 11 ) 
(A. 12) 


The / 2,3 come from the pole lines crossed when shifting the contour. Notice that all the A 
integrals converge rapidly as the T’s decay exponentially in the imaginary direction. 

Let us now take the e —?■ 0 limit of (A.8). Because of the explicit r(e)“^ prefactor, we 
are only interested in obtaining the e poles of the R’s. 


The contour of integration in the double MB integral Ji was chosen serendipitously so 
that it remains 0(1) away from any T function pole as e —)■ 0. Because of this A = 0(1) 
and does not contribute to the leading asymptotics of T. 
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For the second integral we have: 


, ,, r(j)r(-^) r d. ir(-i.)]^r(f+ .) r(| + .-|)r(^ + . + i)r(-j-. + |) 
r(f) J 27ri r(^ + .) r(-|-. + 5) 

(A. 13) 

Here we kept the e dependence only in the last group of factors which become singular as 
e —0 because of the position of the contour ats = —1 + ^ — |. By shifting the contour 
to the right of s = — | + ^ we can arrange so that the new position of the contour remains 
0(1) from any pole. Picking up the contribution of the crossed pole, we get: 


h 


r(!)'r(i)r(-|) 

r(¥) 


r(£/2) + 0(i). 


(A.14) 


Apart from an explicit factor of r(e/ 2 ), in the rest of J 3 we can set e to zero. We get: 


r(f) 


■i) c)r( ^ Ar(^ i.y-rC)PAP Ad)'Fit 


(A. 15) 

Barnes’ lemma was used here to express the integral (the contour in the integral passes 
between s = —d/2 and s = —(i/4, separating the P function poles with positive and negative 
s in the argument). Therefore for Is we get the same answer as for I 2 in (A.14). 


Finally, plugging the found values of 12,3 into (A. 8 ), we obtain (A. 2 ). 


A.2 Gegenbauer polynomial method 


This method expands massless propagators like l/|x — j/p" into powers of |a;|/|?/| or the 
inverse, depending on whether |x| < \y\ or \y\ < |x|, times angular factors which are 
Gegenbauer polynomials. We will use it in a version explained by Kotikov in [79],[80], where 
we refer for notation. The key equation is (3) from [80], which we rewrite as {x = x/\x\)-. 

1 .00 

T = min(|a;|/||/|,||/|/|a;|). (A.16) 

Here the naturally appearing Gegenbauer polynomial is G", while it’s convenient to have 
G;)’s with A = (i/2 — 1 for subsequent integration. To pass from a to A we use (5) of [80]: 


p,, (n - 2fc + A)F(A) V{n + a- fc)F(fc + cr - A) 
k\T{a) F(n + A + 1-A:)F(a-A) ■ 


Substituting back in (A.16) we hnd, in agreement with (3) in [81]: 


(A. 17) 


1 F(A)>^F(n + a) 


|X-|/|2“ (|x||?/|)^ 


r(«)^r(n + A) 


C/^{xy)T^+'^ 2 Fi (n + a, a - A, n + A +1, T'). (A. 18) 
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Finally (and this is Kotikov’s trick) let us represent Gegenbauers as the products of symm- 
metric traceless tensors, eq. (7) of [80]: 


We obtain 


C^{xz) 


2”r(n + A) 
n!r(A) (|a;||^|)"- 


_ _ _ \ ^ ^_ \ ^ pr/« rpn+2k+a 

.y\2a ^ (1x111/1)“+- ^ 

u/« = 2^ 1 r(n + a + fc)(a-A)fc 
n!r(«) fc!(n + A + l), ' 


(A.19) 


(A.20) 

(A.21) 


A.2.1 Integration rules 


In this section parameters a and {3 are arbitrary, not necessarily the same as in our integral 
(A.l). Denote Dx = d'^x/(27r)'^, (n) = pi... /i„, = 6{x‘^—y‘^). We will need the following 

integrals, see (11),(12),(14),(15) of [80]: 
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(n) 


x2"(x — y)'^3 
1 


«... = !/'”> 




^ |^|2(o+/3)-d 
— n I I 
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2 \m / 2\"i+4+fi 
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y ,2 


Here {u = l/( 47 r)'^/^): 

K{n, a, (3) = a;A"'’°(a;, /d), 

BDfi TL\j3 A) 

Li(m,n,a,/?) = a;—-— ’ ’ , , Mi(n,a,/?) = ^ Li(m, n, a,/3), 


m=0 


(A.22) 


(A.23) 


m + a + /S — d/2 
B{jn n\j3 A) 

L 2 (m, n, a, /3) = a;-^, M 2 (n, a, /3) = Y] L 2 (m, n, a, /3), 

m —a + n + a/2 


(A.24) 


m=0 
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while A”’°(q;,/ 3) and -B(m,n|/3, A) are given by 


Al"’°(a,/9) = 


V{d/2 — a + n) 


®n(Q^)no(/3) / _ 

a„(a +/3 — (i/2) ’ ^ r(a) 

r(m + n + /i)(/i - A)^ 


B{m, n, /S, A) = 


m\r{m + n + d/2) 


(A.25) 


It’s possible to perform the summation in the dehnitions of Mi^2, with the result expressible 
in terms of 3-F2’s at z = 1: 


Mi{n, a, (3) 


M2{n,a,(3) 


2^-dT^-d/2r{n + /3) 

~ ri/3)i-2a -2/3 + d)r {lid + 2n)) 

^3 i/j (— d/2 + /3 + 1,n + /i, — d/2 + a + /i; d/2 + n, — d/2 + a + /3 + 1; 1) , 


2i-d^-d/2Y(^n + j3) 

r(/i)(-2a + d + 2n)T {\{d + 2n)) 

X 3-^2 (*i/2 An. — a, — d/2 A /3 A 1, n A /ij d/2 A n, d/2 A n 


(A.26) 


a A 1; 1) • 


A.2.2 Computation 


Let’s start with the integral in y. 


Dy-. 


(A.27) 


We will expand the \x — y\ propagator. (We checked that expanding the \y — z\ propagator 
gives the same result.) After the expansion we get 


All the integrals can be done using Eq. (A.23), producing terms of the form times 

9x:,z or times powers of |x| and \z\. We will then multiply this by 

\x - z\-^'^\z\-^^ (A.29) 

and do / Dz, using (A.22). The arising products can be reduced using (7) of [80]: 

^ (2A)„/[2”(A)„]. (A.30) 

The end result will collapse to an inhnite series of terms all proportional to |x| to the same 
power, and we have to sum the series to get the prefactor. 


Let’s supply the details. Integrating the 9x,y term in (A.28) we get 

/ |^|-2a-2fc-2n|^|-6a+(i+2fc^^^^ 2a - k, 0)9^,z 
n,k ^ 

m 

+ XI n, 2a - k, a)9zA • (A.31) 

m ' 
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Now let’s specialize to a = {d — e)/4. We will only be interested in extracting the leading 
dependence in the e —?• 0 limit. As we will see, the end results is hnite in this limit; 1/e 
poles which appear at the intermediate steps of the computation cancel or get multiplied 
by 0(e) coefficients at later stages. 

It’s interesting to trace how these potential 1/e terms disappear. From (A.28), we can 
expect a 1/e dependence in the integral of the 6x,y term when k = n = 0. In (A.31), this 
divergence is present in the K terms for n = k = t) and in the L 2 term for n = fc = m = 0, 
with the same leading coefficient. Gathering these divergent terms, we have: 


2l-d^-d/2 

r(d/2)e ^ 


—2a I I —2a 


{\z\'^x,z + \x\"dz,x) ■ 


(A.32) 


This grouping makes it easy to see that in the subsequent integration in z these terms will 
be multiplied by an 0(e) coefficient, giving an 0(1) contribution to the hnal answer. This 
can be checked explicitly in the formulas below. 

Analogously, the integral of the 6y^x terms in (A.28) gives: 

J 3a + k + n, a)e,,x 


n.k 
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X 


L2(m, n, 3a + fc + n, a)6z,i 


—QoL-\-d-\-2m\ \—2<y.—2m—2n 


(A.33) 


None of these terms are singular in the e —?• 0 limit. 

Summing (A.31) and (A.33), multiplying by (A.29) and doing J Dz, we obtain:^° 

T = {2TTf<^ fnW^A [K{n, 2a - k, a)M 2 {n, 5a - d/2 - k, a) 


n,k 


+ K{n,3a + k + n, a)Mi{n, 6a — d/2 + fc + n, a)] 

2(m, n, 2a — k, a)Mi{n, 3a + m + n, a) 

— L 2 {m, n,3a + k + n, a)Mi{n, 3a + m + n, a) 

— Li{m, n, 2a — k, a)M 2 {n, 2a — m, a) 

+ Li{m, n,3a + k + n, a)M 2 {n, 2a — m, a)] |. (A.34) 

Among the Mi and M 2 coefficients here, the only one which can become singular as e —)■ 0 
is M2(n, 2a — m,a). The terms with this coefficient came from the integration of the second 
terms in (A.31) and (A.33). The singularity happens only for m = n = 0, and is given by: 


2d_— 


tt-'T (-f+ fc + i)r(f+ fc) 


£Hr(i-|)r(f) [r(f)]qi) 


+ 


k " 


d — Ak d + 2k 


(A.35) 


50r 


The factor (27r)^“ comes from the definition of Dx. 


43 









Mathematica sums the series in A; to a product of T functions, FullSimplify’ed to zero. 

Denote the terms associated with the six lines of (A.34) as i = 1.. .6. We have: 

limT = /i + /2 + t/i + t/2 + >^3 + >^4 5 (A.36) 

e —>0 

T _ ^(1) ^(3) I 

'^k=n=0 ' '^fc=n=m=oh^O ; 

h = y2AC^"^+C^%=m=0U^0, 

■Jl = Y C^^^le^O, 

^(fc,n)7^(0,0) ' 

J2 = y" c'^i.^o, 

•^3 = E, c"*’!"”' 

c(^) + 

fc,n,m:(n,m) 7 ^( 0 , 0 ) 

We omitted as it vanishes for e —?• 0. We also took into account the above discussion 
about the separation of the 1/e poles. The A is hnite for e —?• 0. In I 2 , the individual terms 
are ~ 1/e with coefficients which sum to zero as we mentioned above. So we can drop these 
singular parts and sum the 0(e°) pieces. The A series are manifestly hnite as e —)■ 0. 

We used (A.36) to perform a numerical check of (A.2) for a few values of d. Although the 
representaion does not converge very rapidly, including 200 terms in the single series and 
terms with n, fc, m up to 40 in the double and triple series, we found an agreement within 
one percent or better. This check complements the MB result from the previous section. 




B Conformal invariance of composite operators 

We would like to discuss to what extent conformal invariance of the correlators of (p can be 
used to show conformal invariance of the composite operators, as mentioned in section 4. 

We will use the fact that composite operators can be obtained by fusing (p’s via the 
OPE. The OPE is well established in perturbation theory [24]. Here we will use it in the 
IR where all operators acquire their hxed point dimensions. So we really set A = A*. Let’s 
start by considering a simple example—the operator In the IR region the relevant 

OPE takes the form (A = A^, A2 = 

(p{x)(p{0) = + |a;|^^(p^(0) + ...). (B.l) 

This equation can be used to determine correlation functions of (p^ if correlators of (p are 
known. Indeed consider the correlator 

{p{xi)p{x2)...). (B.2) 

this section all operators will be finite, renormalized, so we drop square brackets from the notation. 
^^The value of the constant coefficient with which appears in this OPE is unimportant for this 
argument, so we set it to one. 
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We then have: 

{^p^{xi)...) = lim \xi 2 \‘^^~^^{{(p{xi)ip{x 2 )...) - \xi 2 \~^^{... (B.3) 

X2^Xl 

This limit should be understood as Xi ^2 becoming much closer to each other than to any 
other insertion. 


Now, Eq. (B.3) implies that correlators of are conformally covariant. Indeed, under 
a conformal transformation x —>■ f{x) both terms in the RHS acquire the same factor 

[J{xi)J{x2)]-^^^^ ^ J{xi )-^^, (B.4) 

where J = \df/dx\^^^ is the a:-dependent scale factor. This is precisely the transformation 
rule of a primary of dimension A2. 


The above argument can be repeated by fusing other pairs of (p’s. This proves that 
correlation functions of an arbitrary number of (p’s and an arbitrary number of (p^’s are 
conformally covariant. 


The next step is to discuss the derivatives of (p^, which appear among the ... terms in 
(B.l). We claim that the coefficients of these terms are such that this part of the OPE 
respects conformal invariance. To prove this consider the three point function 

((p2(l/)(p(a;i)(p(x2)). (B.5) 

By the proven conformal invariance property, it has the well-known functional form [37] 


const. 


\xi — y\^^\x2 — ' 

Expanding this in X 2 — X\, we will get a series of the schematic form 


const. 
xi — 





(B.6) 


(B,7) 


As is well known, matching this series to the one obtained from the OPE, we can £x all the 
coefficients of the 5“(p^ terms, consistently with conformal invariance. Moreover this proves 
that (p^ has zero two point functions with any other operator O appearing in the OPE (p x (p. 
Indeed, a nonzero two point function would give rise to a term scaling as l/|a;i — . 

Assuming as we will, that in a generic interacting theory such as the WE hxed point all 
nontrivial operator degeneracies are lifted, there is no room for such a term in (B.7) except 
if (P is a derivative of (p^. 


Let us then consider the next operator in the (p x (p OPE which is not among the 
already considered operators 1, (p^, and the derivatives of (p^. It will be the operator (p^. 
Its correlation functions can be dehned from those of (p by the equation: 

{ip^{xi)...) = lim |a:i2|^^“^^(((p(a;i)(p(a:2) • • •) 

X2^Xi \ 

-\x,2r^{...) 

- \xi2\~‘^^^^''{P{x2 - xi,d)ip^{xi )...)). (B.8) 
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Here P is the infinite series encapsulating the contributions of and of its descendants in 
the OPE.^^ This part of the OPE has already been proven to be conformally invariant. It 
is then not difficult to convince oneself that all three terms in the RHS of (B.8) transform 
under conformal transformations as expected for a primary of dimension A4. 

It should be clear that the above argument can be continued recursively to higher 
dimensions. At each step we will identify a new operator, show that it’s a primary, and then 
show that contributions of all of its descendants in the OPE are consistent with conformal 
invariance. We then subtract away this part of the OPE, identify the next operator, and go 
on. In the above examples the considered primaries and were scalars, but modifications 
to tensor operators are straightforward. 

The above argument thus determines conformal transformation properties of all opera¬ 
tors in the ip x (p OPE. It is expected that all Z2 even primary operators which are scalars 
or symmetric traceless tensors of even spin appear in this OPE. Z2 odd operators can then 
be accessed via the OPE ip x (p^. Mixed Lorentz symmetry representations can be accessed 
by considering OPEs of operators with spin. This will be somewhat more complicated since 
there are several conformally covariant structures for the three point functions of operators 
with spin. This complication is only technical, and the above argument can be suitably 
extended. We will not discuss this in detail. 


C Relative normalization of the (f> and OPE coeffi¬ 
cients 


In this appendix we investigate the consequence of the nonlocal equation of motion (3.21) 
for OPE coefficients and operator normalizations. We will be working in the IR theory. 
With a small redefinition of the constants involved we can write (3.21) as 

J = C(j)%x) , C = -^{27ry^Wd-aWd+a > 0 . (C.l) 


Let us first investigate the consequences of the equation of motion for three point 
functions. Using (C.l), and Symanzik’s star integral formula [78], we can easily deduce 
that if 


{<^{x)OM02{z)) 




(C.2) 


then 


{,j,^(x)Ody)0,(z)} = 


'1203 


|x — y\ P 


A ,3+Ai-A2 


\X — Z 


A ,3—A1-I-A2 


\y-^\ 


— A^3-|-Ai-|-A2 


(C.3) 


^^Only the first few terms of this series, corresponding to the descendants of dimension < A 4 , will 
contribute in this equation. 
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with Afp3 = d — Afj) and with a relative three point fnnction coefficient given by 
Ai203 _ M 3 _ - A^ + Ai2))T{^{d - A^ - Ai 2))T{A^ - d/2) 

Ai 2* ^ c ’ r(<i-A*)r(i(A* + Ai2))r(i(A*-A,2)) 

Here we introdnced A12 = Ai — A2. Notice that this compntation provides an alternative 
derivation of the scaling dimension A^a in the IR theory, which was derived from the two 
point fnnction in the main text. Notice also that the ganssian limit e —)■ 0 is smooth becanse 
in this limit A is only nonzero if A^ = ±Ai2 — 2k, and the zero in M3 precisely cancels the 
overall 1/go ~ 1/e. 


The above three point fnnctions are rather meaningless withont knowledge of the two 
point fnnctions of the operators involved. The two point fnnction of 0 takes the form 

vO ^ 

The coefficient p is a nontrivial fnnction of e which is hnite as e —?• 0. The two point fnnction 
of 0^ is then given by the analogne of (3.37), which in the cnrrent conventions takes the 
form 

= I rfy I rf. . (C.6) 

Notice that we snbtracted the tree-level part from the two point fnnction of 0, cf. the 
discnssion below eqnation (3.37). The resnlting two point fnnction is therefore proportional 
to p. The integral as written is actnally singnlar, bnt we can view it as a formal expression 
that has a well-dehned and hnite meaning in momentnm space. It then evalnates to 

{^(1)^(0)) = 7010! (C.7) 

X ^ 


with 

^ - A^)r(A^ - rf/2) 

" r(d-A*)r(A*) 


We emphasize that the ganssian limit is again smooth. 


(C.8) 


Snbstitnting A<^ = (d — e)/4 in (C.8) we hnd that M2 < 0 in a hnite interval aronnd 
e = 0, for all physical valnes of d. The coefficient p is then also expected to be negative, so 
that altogether the coefficient of (C.7) comes ont positive. Indeed we expect the LRI to be 
described by a nnitary, or rehection positive in the Enclidean, conformal held theory.We 
see from (3.26) that p is negative to the hrst nontrivial order in pertnrbation theory, which 
is in agreement with expectations. 


Meaningfnl three point fnnctions involve the nnit normalized operators 


0(x) 


1 

a /1 + eV 


0 (x), 



1 

^e^pM2/C^ 



(C.9) 


®^The gaussian theory is unitary as long as is above the scalar unitarity bound d/2 — 1, i.e. if cr < 2. 
This condition is satisfied throughout region 2, see Fig. 1. We expect that perturbing a unitary theory by 
a hermitian operator 0^ with a real coupling will give rise to a unitary theory. 
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in terms of which we hnd that 


•^1203 _ ^71 + 


A 


120 ^/e^pM^ 


(C.IO) 


We see that C drops out, but the relative three point functions still involve p which we can 
only compute perturbatively. 

Notice however that we can take further ratios to get rid of the unknown p. Namely, if 
we consider two OPEs Oi x O 2 and 0[ x C>2, then 


■^ 1203/'^120 
'^1'2'03/-^1' 


2'0 


M' 


(C.ll) 


This equation may be used in a variety of ways in conformal bootstrap analyses involving 
multiple correlators. 

It is interesting to consider the limit cr —?• a* where the LRI_should transition to the 
SRI. The SRI does not have an operator corresponding to our 0^. However from (C.IO) we 
see that there is no decoupling of the (f)^ operator unless e^p = — 1, which would imply that 
0 becomes null and also decouples. Barring this decoupling, we predict a discontinuity in 
the transition from LRI to SRI.^^ 
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